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i-^ ' Abstract 

The LSND and MiniBoone seeming anomalies in neutrino oscillations are usually attributed to 

Ch ' physics beyond the Standard model. It is, however, possible that they may be an artefact of the 

2 . theoretical treatment of particle oscillations that ignores fine points of quantum measurement theory 

^^' relevant to the experiments. In this paper, we construct a rigorous measurement-theoretic framework 

for the description of particle oscillations, employing no assumptions extrinsic to quantum theory. 

The formalism leads to a non-standard oscillation formula; at low energy it predicts an 'anomalous' 

oscillation wavelength, while at high energy it differs from the standard expression by a factor of 2. 

r^ ' The key novelties in the formalism are the treatment of a particle's time of arrival at the detector 

f^ , as a genuine quantum observable, the theoretical precision in the definition of quantum probabilities, 

CO ' and the detailed modeling of the measurement process. The article also contains an extensive critical 

\l^ , review of existing theoretical treatments of particle oscillations, identifying key problems and showing 

lO ' that these are overcome by the proposed formalism. 

o 
o 

1 Introduction 

k>( I In this paper we construct a solid measurement-theoretic framework for the description of particle 
%^ • oscillations [H [2], i.e., a framework that (i) contains no ad-hoc assumptions extrinsic to the general 
^ rules of quantum theory, (ii) allows for a precise correspondence between experimental operations and 

elements of the formalism, and (iii) is sufficiently general to incorporate fine details in the description of 
particle oscillations, when such are needed. The proposed formalism predicts an 'anomalous' oscillation 
wavelength at the low end of the neutrino energy spectrum, and, as such, it might provide an explanation 
for excess of events at low energy in the MiniBooNE experiment and of the LSND anomaly O IH [5]. 

Existing treatments of particle oscillations, as a rule, introduce assumptions extrinsic to quantum 
theory, hence important problems arise. In particular, existing methods contain ambiguities in the rela- 
tion of elements of the formalism to the experimental setup and in their treatment of time and quantum 
probabilities. These problems originate from the fact that particle oscillations touch upon fine points of 
quantum measurement theory (the quantum treatment of the time of arrival, the joint measurement of 
incompatible observables and the role of decoherence and coarse-graining in the measurement process), 
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whose comprehensive treatment had been unavailable at the time the phenomenon was first predicted 

The original intent of this paper was to employ the techniques and ideas from quantum measurement 
theory, in order to resolve the problems above and to place the standard results on particle oscillations 
on a firm conceptual and mathematical footing. Impressively, we find that our method leads to a 
non-standard expression for the wavelength of particle oscillations, i.e., to an expression that disagrees 
with the result that has been obtained by the vast majority of existing approaches to the treatment of 
particle oscillations. We have verified our result to be robust: it is not an artefact of any calculational 
approximation, its derivation involves no ad hoc assumptions and it persists when more details are added 
in the models of the measurement process. The approach we propose is based on a novel treatment of 
the concept of time, as presented in the development of the History Projection Operator theory as a 
new approach to quantum gravity [3 [8]. As far as the theory presented in this paper, the distinction 
between the kinematics and dynamics on the definition of time resulted to a new treatment of the time 
of arrival as a genuine quantum observable [9] . 

Claims of non-standard oscillation formulas have appeared sporadically in the literature — for exam- 
ple, Refs. [ini H] — and they have been repeatedly refuted in favor of the standard result [I2l [I3l [HI [I5l 
116] . We do agree with the critique of such non-standard formulas, in the sense that the methodologies 
through which they are produced can be found wanting. However, we claim that the standard deriva- 
tions are also wanting on the same grounds. To show this, we have undertaken an extensive critical 
review of existing theoretical treatments of particle oscillations (Sec. 4), identifying their shortcomings 
and showing that these shortcomings are, indeed, overcome with the method we develop here (Sec. 5.3). 

One argument in favor of the non-standard result obtained here is that its derivation is not plagued 
by any mathematical ambiguities or conceptual inconsistencies. More importantly, our expression for the 
oscillation wavelength can be experimentally distinguished from the standard one, at least in principle. 
For neutrinos, the standard expression for the oscillation wavelength A at energy E is 

where Am'^ is the difference of squared neutrino masses. The non-standard oscillation formula that has 
been obtained before differ from Eq. ([I|) by a factor of two, i.e., 

A - ^^^ (2) 

In absence of an independent precise measurement of neutrino masses Eqs. ([1]) and ([2]) cannot be 
experimentally distinguished. In contrast, here we find 

Am2(l-|i)' ^^ 

where eth is the energy threshold for the process through which the neutrino is detected. While for 
E » eth, Eq. ([3]) coincides with Eq. ([2]), for £■ near eth, Eq. ([3]) predicts that the ratio X/E is not a 
constant but increases with energy. 

Hence, the method we develop here leads to the following predictions. 

• At low energies, the oscillation wavelength depends strongly on the value of the threshold energy 
eth, hence, it may vary according to the process through which neutrinos are detecteco. In other 
words, data subsets for different reaction channels may exhibit different oscillation patterns. 



^To be precise, Eq. ((3| is obtained from the study of reactions in which the neutrino is annihilated, i.e., charged current 
interactions — see. Sec. 5.2. Neutrino detection through scattering processes can also be studied in the present framework 
and it will be treated in a different publication. 



• For a given reaction channel, the expected number of events as a function of energy will differ 
strongly from what is predicted by the standard oscillation formula at low energies. This difference 
can be seen heuristically from the plots in Fig. 2 (Sec. 5.3). 

The prediction of 'anomalous' neutrino oscillations at low energies might provide an explanation for 
the excess of neutrino detection events in the LSND [3j and MiniBooNE |4j experiments. In particular, 
our approach makes no recourse to physics beyond the Standard Model and introduces no additional 
undetermined parameters. Hence, it can be directly tested against the standard expression for particle 
oscillations in terms of better fit to existing data. 

We will refer to the formalism developed here as the quantum measurement approach to particle 
oscillations. Its main points are the following. 

• Particle oscillation experiments do not correspond to measurements taking place at a predetermined 
moment of time. The time of arrival t of the particle at the detector is a genuine physical observable 
that should be described in terms of quantum theory. In particular, the time of arrival is to be 
distinguished, both physically and mathematically, from the time parameter t of Schrodinger's 
equation [7]. 

• For a quantum measurement the most general consistent assignment of probabilities is through the 
use of Positive-Operator- Valued-Measures (POVMs) — see Sec. 3.1.3 for definition. The probability 
densities for particle oscillation are represented as Pa{t, L) = Tr (poIla{t, L)] , in terms of a family 

of POVMs tla{t, L), where t is the time of arrival, defined in terms of the detection of a particle 
produced through a channel of flavor a, and L is the distance from the production region. Hence, 
we describe particle oscillation experiments as joint measurements of the time of arrival t and 
the position L. This is a key difference of our approach from existing ones, and it leads to a 
non-standard expression for the oscillation wavelength. 

• The family I[a{t,L) is constructed in terms of the detection of product particles, i.e., in terms of 
directly observable quantities. Moreover, the method for constructing Ilait,L) involves only the 
Hamiltonian of the system and the specification of the measured variables. In particular, it does 
not depend on properties of the initial state and contains no assumptions extrinsic to quantum 
theory. 

The paper is structured as follows. Sec. 2 describes our motivation in undertaking this work, 
in order to correct certain conceptual and mathematical inadequacies in existing treatments of particle 
oscillations. We identify theoretical assumptions that enter into existing derivations that are not justified 
by the rules of quantum theory and, sometimes, they are erroneous. In that section we also present the 
main ideas of our approach. The presentation is largely non-technical, and it requires little background 
to the issue. 

Sec. 3 provides background material necessary for the development of our approach. In particular, 
we focus on quantum measurement theory providing a review of main results that are relevant to our 
critique of existing approaches to particle oscillations. 

Sec. 4 presents a specific and more technical critique of existing approaches to particle oscillations, 
in light of the inaterial presented in Sec. 3. 

The development of a new framework for particle oscillations is contained in Sec. 5. In Sec. 5.1 
we construct the necessary tools in the context of an abstract Hilbert space formalism. In particular, 
we develop a methodology for defining probabilities for general measurement schemes, where the time 
of detection is treated as a physical observable. Particle oscillation experiments constitute a special 



case of such schemes, and they are treated in Sec. 5.2, where we derive the formula for the oscihation 
wavelength. In Sec. 5.3 we undertake a critical analysis of our method, demonstrating that it resolves the 
problems we identified in existing approaches to the issue. In Sec. 5.4 we further develop the formalism, 
applying it to more general settings, and we demonstrate that our results are robust. 

2 Motivation for a quantum measurement approach to particle oscil- 
lations 

The problems in the existing treatments of particle oscillations arise from their identification of prob- 
abilities as the modulus squared of matrix elements of the evolution operator e~*^*. While this is the 
most common probability assignment in quantum theory, it is not universal. In particular, it does not 
apply to particle oscillations: the resulting probabilities depend on the time t and there is no rigorous 
way to relate them to what is actually being measured in an experiment. When particle oscillations were 
first predicted, the researchers did not have the benefit of a mature quantum theory of measurement 
that could provide a solution in terms of first principles. As a result, the predictions relied on heuristic 
arguments and assumptions that are not warranted by quantum theory. These created ambiguities in 
the relation of the formalism to experimental procedures. 

Later work on the quantum theory of measurement provided a precise language for the treatment 
of the subtle issues that enter the theoretical description of particle oscillations. As it turns out, the 
more precise treatment yields information inaccessible from the heuristic ones, namely, the prediction 
of 'anomalous' oscillation wavelength at low energies that could be relevant for the explanation for the 
LSND and MiniBooNE experiments. 

In section 2.1, we provide an extensive discussion of problems in existing treatments of particle 
oscillations, explaining why they arise and how they affect the existing treatments. In section 2.2, we 
describe the main ideas of our approach, and how they resolve these problems. 

2.1 Main problems in particle-oscillations theory 

We proceed to an analysis of the following issues: (i) the improper treatment of the time of arrival, (ii) 
the lack of a first principles derivation of the relation between amplitudes and probabilities, and (iii) 
the improper elimination of the time variable in the probability assignment. In our opinion, these issues 
constitute the most important problems in existing treatments of particle oscillations. 

2.1.1 The time-of-arrival issue 

Current methods for the treatment of particle oscillations do not establish a sharp correspondence 
between mathematical formalism and experimental operations. The reason is that these methods conflate 
two distinct notions of time as it appears in quantum theory: time as a parameter of evolution in 
Schrodinger equation and the time of arrival (or time of event, or time of transition) [7]. 

Time as a parameter of evolution, like in Schrodinger equation, is an essentially classical parameter 
[T7j . i.e., it does not correspond to a Hilbert space operator; hence, it is not an 'observable'. On the 
other hand, the time of arrival is a variable defined in terms of coincidences of events: we observe an 
event (e.g., the creation of a lepton in a detector) and we correlate its occurrence with the reading of an 
external clock. The time of arrival is, therefore, a physical observable. 

In existing treatments of particle oscillation the distinction above is ignored. This results into the 
definition of probabilities in terms of the Schrodinger time parameter which, unlike time of arrival, is not 
a physical observable. To obtain physical predictions from these probabilities it is necessary to remove 



this time dependence; this is achieved by procedures that are not justified from first principles — see Sec. 
2.1.3. 

The time of arrival is defined in terms of the physical event that a particle is detected by a measuring 
device. Since particle detection is a quantum process, time of arrival is a genuine quantum observable 
associated to physical processes on the measured system. (See, Sec. 3.3 for details and references.) As 
such, time of arrival may exhibit interferences. In a hypothetical experiment where a source prepares 
particles in a superposition of states with different mean momentum, the times of arrival recorded at 
a detector, at distance L from the source, will be distributed according to an oscillatory pattern (i.e., 
interference) around a mean value of arrival time. 

2.1.2 The physical interpretation of amplitudes 

The vast majority of the theoretical studies of particle oscillations concentrate on the evaluation of 
amplitudes between certain vector states on a Hilbert space that describes the physical system. For 
example, these states may be 'eigenstates' for flavors a and /3, leading to amplitudes of the form 

AaMt^-L) = (a,x = 0|e-*^*|/3,x = L). (4) 

In Eq. dH) the states |x, a) refer to particles localized around the point x and they correspond to detection 
of particles characterized by flavor a. In other approaches that employ quantum held theory, analogous 
amplitudes are computed between states of incoming and product particles, rather than neutrinos or 
neutral bosons, as in Eq. ^. The dependence of the amplitude on the flavor indices arises through the 
flavor dependence of an interaction Hamiltonian. The critique in this section applies to those approaches, 
as well. The critique refers to a structural point of the quantum formalism (i.e., the relation of matrix 
elements of e~ to the probabilities of measurement outcomes), and not to properties of specific physical 
systems or calculational techniques. It also applies to cases that the amplitude is defined in terms of 
elements of the 5-matrix rather than matrix elements of the finite-time propagator e~*^*. 

Having computed an amplitude of the form Eq. (j4j), the modulus square \Aa^p{t,'K)\'^ is then inter- 
preted as a measure of probability, from which physical predictions about particle oscillation experiments 
are extracted. In Sec. 4, we will examine the various methods employed to extract physical predictions 
from such amplitudes. 

A key question then arises: in what sense are amplitudes of the form Eq. (j4]) relevant to the description 
of particle oscillation experiments? In particular, Eq. ^ defines the matrix elements of the unitary 
evolution operator between states of the form |a,x). How does this relate to probabilities relevant to 
particle oscillation experiments? 

The squared amplitude \Aa^p{t,'L)\^ defines probabilities for the measurement of a particle's posi- 
tion, at a specific, predetermined and sharply defined moment of time t. Besides the idealization involved 
in talking about measurements taking place at a sharply defined moment of time, the main problem in 
the use of the squared amplitude \Aa^p{t,'L)\^ is that in particle oscillation experiments the timing of 
the detection event is not predetermined by the experimentalist. What is controlled is the location of 
the detector. In other words, particle oscillation experiments do not fall in the domain of validity of the 
expression |^fj_!.^(t, L)p for probabilities of measurement outcomes. We shall elaborate on this point in 
Sec. 3.1. 

In order to construct a proper probability assignment that fully reflects the setup of particle oscillation 
experiment, it is necessary to work within the broader context of quantum measurement theory [18 1 119 ^ 
I20j . The latter allows for probability assignments more general than squared amplitudes (for example, 
probabilities for joint measurements [18], continuous-time measurements [21], weak measurements [22], 
time-of-arrival measurements [231 124^ 19]. and so on). In particular, it is necessary to derive the relevant 



probability from first principles, taking into account the physics of the detection scheme; we shall provide 
such a derivation in Sec. 5.1. 

In most physical applications, a first-principles derivation of the probability assignment through 
quantum measurement theory makes little difference to the theory's predictions, hence, it might be 
dismissed as too fine to be relevant. However, this is not the case for particle oscillation experiments, 
in which the measured quantities are directly correlated with tiny microscopic variations (oscillations) 
in the time of arrival of the particles. A proper treatment of the probabilities does affect the physical 
predictions of the theory. 

2.1.3 Elimination of the time variable 

The problems in the use of the squared amplitude \Aa-^i3{t, L)p for the description of particle oscillations 
have been noticed early on in the theory, as the presence of the time t in addition to the position L 
appears superfluous when relating the theoretical expression to experimental results. However, this issue 
has not led to a reconsideration of the starting point, that is, to the derivation of the probabilities for 
particle oscillations from flrst principles. Instead, the relevance of the squared amplitudes has been taken 
for granted, and various procedures have been invoked in order to eliminate time from the probabilities. 
Clearly, since the starting point is accepted without justification, any remedy to the consequent problems 
can be introduced only ad hoc, i.e., without justification on the basic principles of quantum theory. 

Here we describe the two most common such procedures, for eliminating time from squared ampli- 
tudes. 

Use of the classical arrival time. The oldest approach consists in using |^Q,^^(t, L)p as a proba- 
bility density for the detection, by substituting t with a classical time of arrival t^i = L/v; L = |L| is the 
distance from the detector to the source and v the mean velocity of the wave-packet. 

The obvious problem with this method is that it is mathematically inconsistent: the time t, that a 
measurement takes place, is substituted by an 'average' time-of-arrival, which depends on the system's 
initial state, through v. The resulting probabilities are then not linear functionals of the initial state 
of the system p, as required by quantum theory. Hence, this solution can only be a heuristic device, 
perhaps valid for a specific class of initial states. 

The method above could be viewed as a kind of semi-classical approximation. As such, it can only 
be justified by demonstrating that it is a good approximation, i.e., by comparing with a full quantum 
mechanical treatment of the problem. Heuristic arguments based on estimations of uncertainties have 
pitfalls, and they are no substitute for a full calculation. 

Integration of probabilities over time. A second commonly employed procedure for the elimination 
of time t from the squared amplitude |^a-5>/3(^)L)p is integration. The related argument is that, since 
time is not a directly measured quantity in the experiment, we integrate over the single-time probabilities 
\Aa^li{t,'L)\^. We then obtain a probability density for position alone, i.e., 

p^^pi^L)^ [ dt\A^^p{t,\.)\^, (5) 

Jo 

where T is an integration time for the experiment. 

The argument above is invalid for the following reasons. First, the squared amplitudes \Aa-¥is{t, L)p 
are probability densities with respect to position and not with respect to time, and thus time-integration 
does not lead to meaningful probabilities. This can be seen even on dimensional grounds: a density with 
respect to time as well as space ought to have dimensions \time\~^ x \length\^^ , while the squared 



modulus |^Q,_>^(i, L)p has dimensions [length]~^ . Moreover, it does not transform as a density under a 
time-rescaling t — )• at. 

This obvious problem is sidestepped by a procedure invoked for the solution of another problem, 
namely, the fact that /q dt|^cf_>^(t,L)p is not normalized to unity. One is then forced to normalize 
it by dividing by a constant C, for example, C = ^2/3 I d^^ lo l-^a-s>/3(*)L)p. However, C depends on 
the initial state po, hence, the probability assignment is not a linear functional of po, thus violating a 
fundamental property of quantum mechanics. For an elaboration of these points, see. Sec. 4.2.3. 

The mathematical and conceptual consistency of this procedure apart, there are also physical prob- 
lems with the integration over time in Eq. ([5]). Suppose we were able to properly define amplitudes 
Ba^p{t, Jj), different from Eq. Q, in which the parameter t stands for the time of arrival which is a 
genuine observable (we shall see how this is done in Sec. 5.1). Would we be justified then in defining a 
probability density for position 

Pa^l^)- I <it|ea^/3(t,L)|2, (6) 

if we assume that the time of arrival is not directly observable? In classical probability, the answer would 
be affirmative. However, the defining feature of quantum probability, clearly manifested in the two-slit 
experiment, is that alternatives are summed over at the level of amplitudes. Therefore interference 
terms appear. Only if a macroscopic distinction of alternatives is possible are the interference terms 
suppressed. Hence, Eq. ^ for the probability would be justified only if time of arrival is measured very 
finely. In particular, Eq. ([6]) is justified only if the resolution of the time measurements captures the full 
dependence of Ba^p{t, L) on time. Unless there is no time measurement — -or if the time measurements 
are too coarse to distinguish oscillating terms in the distribution of the time of arrival — then one should 
expect that the relevant probabilities are of the form 

p„_,;3(L)~| / dtfi,^;3(t,L)|2, (7) 

Jo 
where interference terms will be in general present. 

In other words, even if all the mathematical problems associated to integration over time at the level 
of squared amplitudes were resolved, Eq. ([6]) does not follow from the basic principles of quantum theory, 
unless one postulates a decoherence mechanism that destroys all interferences for the time of arrival in 
the amplitudes Ba^p{t.,'L). 

2.2 The main ideas of the quantum measurement approach to particle oscillations 

In Sec. 2.1, we argued that the existing methods for particle oscillations suffer from problems of consis- 
tency with the basic principles of quantum theory. These problems can be concisely summarized in the 
following statements. 

1. A key issue in the theoretical study of particle oscillations is the definition of the time of arrival 
as a genuine quantum variable. In existing treatments, time of arrival is either conflated with the 
time-parameter of Schrodinger's equation, or it is treated classically. The former case contains 
a conceptual error that leads to a discrepancy between the formalism and the experiment. The 
latter case involves an approximation, whose validity can only be ascertained by comparison with 
the results of a full quantum description. 

2. There is no first-principles derivation of the probabilities that pertain to particle oscillation ex- 
periments. The modulus square of the amplitude Eq. (j4]) is introduced without justification as a 
probability measure relevant to particle oscillation experiments. Ad-hoc procedures are employed 
in order to extract physical predictions. 

7 



These issues are resolved in the present paper: the quantum measurement approach to particle 
oscillations is based on two key ideas which can be concisely summarized as follows. 

1. The distinction between the time of arrival and the time parameter of Schrodinger's equation is 
the groundwork of our approach. The importance of this distinction originates from previous work 
on the concept of time in quantum theory [71[8], as we stated in the Introduction. The formalism 
presented here strongly separates between those two different notions of time, by representing 
them by different mathematical objects. In particular, the formalism treats the time of arrival 
as a physical observable, directly relevant to particle oscillation experiments and, it provides a 
probabilistically consistent treatment of its quantum fluctuations. 

2. Treating the time of arrival as a quantum observable, implies that the probabilities relevant to 
particle oscillations correspond to the joint measurement of two observables: the time of arrival t 
and the location of a detection event L. In quantum theory, probabilities for joint measurements 
cannot, in general, be obtained by squared amplitudes. However, they are described, with full 
consistency and rigor, in terms of Positive Operator Valued Measures (POVMs) |18t I20j. Hence, 
we obtain the probabilities relevant to particle oscillation experiments through the construction of 
a POVM nQ,(t,L) that describes the probabilities for the detection of a-flavor events at location 
L and time t. The corresponding probabilities are then of the form 

Pa{I^) = Trn [ma{t,L)] , (8) 

where po is the initial state of the system. Since in Eq. ([8]) the time of arrival t is an observable, 
it can be consistently integrated out, leading to a probability density that is a function of L alone. 
The oscillation wavelength Eq. ([3]) is then obtained. 

The implementation of the ideas above involves the construction of amplitudes Aa{t, L) that depend 
on the time of arrival rather than the time parameter of Schrodinger's equation. These amplitudes 
are then employed for the construction of the POVM IIa(t,L) using a method that was developed in 
Ref. [9]. In fact, we do not construct a single POVM for the task at hand but we provide a general 
methodology for their construction. The results, therefore, do not depend on simple idealizations; there 
exists a systematic procedure that allows us to incorporate as many details in their derivation as needed. 

The measurement process is described in terms of the detection of particles produced by reactions 
of the oscillating particles at the detector, hence the use of quantum field theory is necessary. We also 
take into account the fact that the detector is a macroscopic system with a large number of degrees of 
freedom. In spite of the complexity involved in a realistic description of the measurement process, our 
main result Eq. ([3]) for the oscillation wavelength is robust, and it does not depend on variations in the 
modeling of the measurement. Moreover, it does not require any fine-tuning in the properties of the 
initial state. 

3 Background 

This section consists of two parts that address issues important for the critique of existing approaches 
to particle oscillation. It also serves as a theoretical background for the formalism we develop in Sec. 
5. In particular, in Sec. 3.1 we review the main ideas of quantum measurement theory, placing special 
emphasis on measurements of von Neumann type, explaining why they are not suitable for the description 
of particle oscillations and discussing the properties and wide applicability of POVMs. In Sec. 3.2 we 
briefly review the quantum treatment of time of arrival in the language of POVMs, showing that, like 
any other quantum observable, the time of arrival may be characterized by interferences. 



3.1 Quantum measurement theory 

3.1.1 von Neumann measurements 

First, we study the simple model about the measurement of a variable corresponding to the self-adjoint 
operator A, with generalized (continuous) eigenvalues a. This class of models originates with von Neu- 
mann |251 126) — see also |20j — the variation we present here is from Ref. |27J. 

To this end, let Tig be the Hilbert space of the microscopic system and Tiapp the Hilbert space 
corresponding to the degrees of freedom of a macroscopic apparatus. The values of the observable A 
of the quantum system are correlated to the values of an operator X on Tiapp- X is usually called the 
'pointer variable', its values are accessible to macroscopic observation. Hence, its spectrum is highly 
degenerate. We denote the generalized eigenstates of X as \X,X); A labels orthonormal bases in the 
degeneracy subspaces of X. 

Hg is the Hamiltonian operator describing the self-dynamics of the quantum system. The interaction 
between the measured system and the apparatus is described by a Hamiltonian of the form Hint = 
f{t)A(SiK. K = X^A J dkk \k, X){k, A| is the 'conjugate momentum' of X, i.e., the generator of translations 
for the pointer variable X, where {X,X\k,X') = -h=e^^'^Sx\'. The function f(t) in Hint describes the 
'switching-on' of the interaction between system and measurement device. 

For simplicity, we assume that the self-dynamics of the apparatus during the interaction are negligible, 
so that the total Hamiltonian is Htot = Hg + Hint- The corresponding evolution operator is U{t) = 

rt 

Te~^ JO ^ tot(s) ^ where T denotes the time-ordered exponential. 

The microscopic system is prepared in a state pQ and the initial state of the apparatus is |^o)- The 
initial state of the combined system is then po x |^o)(^o|! and it evolves into ptot{t) = U{t)poU^t), 
under the action of the unitary operator U{t). 

With the above assumptions, the probability distribution pt{X) over the value of the pointer variable 
at time t is 

MX) = Y^iX, X\TrHjtotit)\X, X) = Tths [po^X)] , (9) 

A 

where n(^) = Yl,\ C{X, X)C'^{X, X) is a family of positive operators on Hs, and where 



flk 

C{X,X)= I ^^e''^''{k,X\^o) 



(10) 

A common approximation is to assume that the function /(t) is narrowly concentrated around a 
value t = i, i.e., f{t) ~ 5(t — f). In that case, 

PtiX) = Tths {e^'^'poe'^'wiX - A)) , (11) 

where 

w{X)=J2\MX,X)\^ (12) 

A 

and w{X — A) stands for the positive operator J da w{x — a)\a){a\. The function w{X) determines the 
correlation between the pointer variable X and the values of the operator A. 

Eq. (jlip corresponds to the probability for a measurement of the observable A, at a single moment 
of time i [28]. While the degrees of freedom of the apparatus have been included in the description, the 
end result for the probabilities is expressed solely in terms of operators on the Hilbert space Hs of the 
microscopic system. 

9 



3.1.2 Coarse-graining and measurement uncertainties 

Next we elaborate on the physical origin of uncertainties and the role of coarse-graining in quantum 
measurements. These quantum measurement theory issues are crucial for the theoretical treatment of 
particle oscillations. 

The accuracy by which the pointer variable X determines the value of A is given by the weight 
function w{X) in Eq. (J12p : w{X) depends on the properties of the measurement device. The resolution 
in the measurement of A is at least of order 5X, where 5X is the mean deviation of X for the state 
l^o)- In general, 5X is of macroscopic magnitude, since the apparatus has a large number of degrees of 
freedom. For example, 8X cannot be smaller than the thermal fluctuations of X. 

A measurement is a macroscopic amplification of a microscopic event. The inability to distinguish 
between values of ^ at a scale smaller than 6X is an objective fact of the system. There is no way in 
standard quantum theory — without introducing hidden variables — to infer properties of a microscopic 
system at a scale smaller than the resolution 5X |291 [30l [31} I32j . In quantum theory, microscopic 
system and macroscopic apparatus become entangled during measurement and the association of definite 
properties requires a sufficient degree of coarse-graining. 

Eq. (jlip reduces to the usual probability assignment, in terms of squared amplitudes only when 
sufficiently coarse alternatives are considered^ i.e., at scales much larger than 5X. In that case, the 
weight w{X — A), in Eq. (fTT]l . can be substituted by a delta function and Eq. (fTTI) corresponds to a 
squared amplitude |(V'o|e~*^*|a)p. 

Measurements of von Neumann type are not universally applicable because they involve special 
assumptions and idealizations. In particular, the description of time in von Neumann measurements is 
artificial: we assumed an interaction between microscopic system and apparatus determined by a specific 
profile function f{t). The timing of the interaction between microscopic system and apparatus is not 
controlled by the experimentalist and it is not independent of the initial state. For example, the time of 
arrival of particle at a detector, located at distance L from the particle source is not predetermined by 
the experimentalist and it varies according to the energy distribution of the incoming particles. In this 
setup, the parameter controlled by the experimentalist is the distance L. Hence, an accurate modeling 
of such measurements should involve a profile function for position rather than time. 

The idealization in the treatment of time in von Neumann measurements is a good approximation 
when the precise timing of the event is not significant to the measurement outcomes. However, for 
measurements where the time of detection is a physical variable (particle oscillations), a description in 
terms of von Neumann measurements is insufficient. 

A realistic measurement scheme has a finite duration that corresponds to the time it takes for the 
pointer variable to settle to a definite value, correlated to a property of the microscopic system. This 
time-scale is determined by the physics of the apparatus, it typically lies in the macroscopic domain, 
and it sets a lower limit to the temporal resolution of the apparatus. 

In order to provide an estimation for the duration of the measurement we generalize the model 
presented in Sec. 3.1.1, by including a Hamiltonian term Happ for the apparatus's degrees of freedom, so 
that the total Hamiltonian is H = Hs<S'l + l0Happ+Hint- A good pointer variable X must be stable with 
respect to the self-dynamics of the apparatus, so that evolution under Happ does not significantly affect 

the recorded values after the measurement's completion. This means that e*^"''^*Xe~*''^''pp*|^o) — -^l^o)) 
or more precisely 

^(f) := \{yi,Q\e'^apptj^^-iHappt _ x|^o)| « |(^o|^|*o)| (13) 
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and a similar condition must hold for the operator K that generates translations in X. The condi- 
tion above can be satisfied if X and K are coarse-grained quasi-classical variables, i.e., if they have 
large degeneracy eigenspaces that remain approximately invariant under time evolution. In general, the 
evolution of quasi-classical variables involves noise [32l [331 EH EH ESI EZ], so ^(t) in Eq. (fT3]) can be 
modeled by a classical stochastic process. When the effects of the noise can be ignored in the evolution 
of the system, the previous analysis passes unchanged and Eq. (Ilip is obtained. However, this involves a 
specific scale of coarse-graining both for the uncertainties 6X and for the duration r of the measurement. 

The effect of noise, such as ^{t) above, in the evolution on the quasi-classical variables is often 
studied within the quantum theory of open systems [181 EH]- The ignored degrees of freedom of the 
apparatus are then treated as an environment and they are modeled by a thermal reservoir — see, Refs. 
|29^ [39l HOl W\\ I42j . It is natural to expect that the characteristic timescale r of such measurement 
corresponds to the time that it takes a variable X to settle to an asymptotic value, characterized by 
purely thermal fluctuations. Therefore, r is of the order of the relaxation time Trei of the reservoir. 

The limit above is essentially classical, in the sense that it arises from the statistical fluctuations of 
a system with a large number of degrees of freedom. There is a more stringent, unavoidable lower limit 
for the resolution time scale r arising from the quantum nature of the pointer variable. This limit is 
the localization time tioc, i-e., the time it takes the pointer variable to establish classical behavior due 
to its interaction with the remaining degrees of freedom [361 1^ [Bj • Modeling the pointer variable by 
an oscillator undergoing quantum Brownian motion in a thermal reservoir of temperature T [38^ B5l B5] , 



one finds a localization time of order w— ^, where T^ei is the relaxation time [UlIlllIlT]. Such systems 
are characterized by a time scale separation T^^ « tioc « T^ei- For T = 300"-?^, this implies that 
'Tioc >> lO^^^s. We note that the temporal resolution in a solid state detector is much coarser than 
these lower bounds: it is at best of the order of 10~^s. 

Models such as the above are simplistic in comparison to the physics of actual particle detectors. 
However, they serve to identify the irreducibly quantum origin of the temporal resolution in a measure- 
ment scheme, namely, that the pointer variables can be assigned definite classical properties only if a 
substantial degree of coarse-graining is involved. 

3.1.3 Positive-Operator- Valued Measures 

In Sec. 3.1.1, we constructed the positive operators II(-'^) for a measurement device by considering 
a single pointer variable X, whose values are correlated with one quantum variable A. In general, a 
measurement may involve many pointer variables X* each correlated with a different property A^ of the 
quantum system. Again, the probabilities for the measurement outcomes can be expressed in terms of 
a POVM. Here, we briefiy describe the basic properties and uses of POVMs in quantum measurement 
theory. 

Let % be the Hilbert space of a quantum system and let a* be the values of observables correlated to 
pointer variables of a measurement device. The possible outcomes of the measurements span a space Vt. 
A POVM is a map that assigns a positive operator Piu on % to each (measurable) subset U of fi, such 
that (i) HuiVJU2 — n(7i +n(/2, for Ui,U2 C ri, and (ii) Hq = 1. Usually, the POVM H^ can be expressed 
as an integral over U in terms of operators n(a), a G il, such that Hjj = J^ daH{a). 

The most general probability assignment for measurement outcomes in Q , compatible with the rules 
of quantum mechanics is 

p{U) = Tr{poU.u), (14) 

in terms of a density matrix pQ that incorporates the initial state and preparation of the system and 
some POVM H^/. If a procedure provides probabilities for a measurement outcome that cannot be 
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brought into the form Eq. (|14p . then these probabihties are incompatible with quantum theory. They 
may approximate a vahd quantum mechanical expression (for a special class of initial states), but they 
are not themselves fundamental and cannot be taken as a standard for the theoretical description. 

3.2 The quantum description of the time of arrival 

This section contains a brief review on the quantum treatment of the time of arrival, with an emphasis 
on its description in terms of POVMs. The key point relevant to particle oscillations is to demonstrate 
that the probabilities associated to the time of arrival, like those of any other quantum observable, may 
exhibit quantum interferences. 

3.2.1 Key issues in the quantum description of the time of arrival 

An idealized description of a time-of-arrival measurement is the following. At time t = a source emits 
particles prepared at an initial state po- At a macroscopic distance L from the source a particle detector 
is located. An external clock determines the time t coincident with each detection event. Given a 
large number of events, a probability distribution for the time of arrival is constructed. The issue is to 
determine this probability distribution from first principles in quantum theory. 

The difficulty in describing the time of arrival in a quantum system stems from the fact that given 
a solution to Schrodinger's equation ip{x,t), its modulus square \^p{x,t)\'^ is a probability density with 
respect to x at each time t, and not a probability density with respect to t. Moreover, there is no 
operator representing time in the system's Hilbert space. In that case one would use the Born rule to 
determine a probability density for the time of arrival. The existence of a time operator T, conjugate 
to the Hamiltonian (so that the Hamiltonian generates time-translations: e Te" = T + s), is ruled 
out by the requirement that the Hamiltonian is bounded from below [481 I49j . For some systems, one 
may still define a quantum time variable, by choosing some degrees of freedom of the system as defining 
an internal 'clock'. However, clock times fail to be conjugate to the Hamiltonian of the system, with 
the result that they do not forward under Hamiltonian evolution. In other words, quantum fiuctuations 
invariably force clocks to move occasionally 'backwards in time'. 

There are various approaches for the description of the time of arrival |241 [50l [52] , and related prob- 
lems, e.g., transversal time in tunneling [53]. Some of the approaches go beyond the standard quantum 
mechanical formalism by employing properties of 'paths' identified either through path integrals, or the 
Wigner function, or Bohmian mechanics. In this paper we are interested in well-defined probabilities for 
the time of arrival in concrete measurement situations, which means the identification of a POVM tl{t) 
for the time of arrival O [23l Ell ETJ [58]. We also note that the method we adopt in this paper for the 
construction of a POVM for the time of arrival shares many common features with the studies based on 
the decoherent histories approach pU \55[ [56] . 

3.2.2 Time-of-arrival interferences 

Next, we present Kijowski's POVM [23j for the time of arrival of a non-relativistic particle of mass M in 
one spatial dimension. Given an initial state iV'o) at time t = 0, Kijowski's POVM assigns the following 
probabilities for the time of arrival of the particle at a detector located at x = 

"« = I r * i^L) '" ^-""""■■^■iP)\' + I /I * {^) "^ e-''/-'«P)P. (15) 

The POVM Eq. (jlSp is uniquely constructed by the requirements of Galilean invariance, of invariance 
under parity and time-reversal transformations, and of having the correct classical limit. It is defined 
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on wave functions with support on either positive or negative values of momentum, and the values of t 
extend to the whole real axis. It can be extended to all states if we also include the alternative J\f of no 
detection [9], as p{J\f) = 1 - /^p(t). 

In order to demonstrate that the POVM Eq. (|15p may lead to interferences in the time of arrival, we 
evaluate Eq. (jl5[) for an initial superposition of two Gaussian wave packets with the same mean position 
—L, but different mean momentum pi and p2 



Mp)=\^ 



a-^'^' 



g-^(p-Pi)2+ipL _^_ ^-^{p-P2f+ipL 



(16) 



We assume that a\pi — P2I >> 1 so that the two wave packets do not significantly overlap. 

For api » 1, ap2 » 1, and ignoring dispersion, the leading contribution to the probability 
distribution obtained by the POVM Eq. (jlSp is 



1 

Pit)- 



2TTMa 
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+VP1P2 e ^^^ '^^ ' ' ^ '^^ ' ' ^ cos (pi -P2){L 1) 
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(17) 



i.e., there are oscillations of frequency to = 2mI^i ~ -Psl extending for a time interval with width of 
order At = LM\pi — P2\/{Pip2) around the mean value of arrival time td = ^^M ■ If the momentum 
difference Ap = \pi — P2I is much smaller than the mean momentum p = ^{pi + P2)) then io = pAp/M 
and At = MLAp/p^ . Clearly At/td = Apjp << 1, but this has nothing to do with the measurability 
of the interferences in the time of arrival. The interference are observable if At >> r, where r is the 
temporal resolution of the detector. For measurements of the time of arrival, see Ref. ^J and references 
therein. 

The discussion above demonstrates that there exist regimes in which the semiclassical description of 
the time of arrival is inadequate. A common theme in the study of particle oscillations is to identify the 
relevant probabilities by substituting into the square amplitude Eq. (j4]) the classical value of time of 
arrival at the detector. The meaning of such an expression apart, this procedure may gravely misrepresent 
the probabilities of the system. Indeed, for the initial superposition state Eq. (fTBj) substituting, for 
example, t = id into the square amplitude |?/'(t,a;)p, leads invariably to an expression that does not 
exhibit an interference pattern. 

4 Critique of existing theoretical treatments of particle oscillations 

In this section, we undertake a more detailed critique of existing approaches to particle oscillations. 
First, we present the existing approaches for the treatment of particle oscillations, separated into three 
broad categories as in the reviews |59l [6OI [611 [62] , namely, plane wave methods, wave packet methods 
and quantum field theoretic methods. We then proceed to a critique of their assumptions and formalism, 
elaborating on the issues sketched in Sec. 2. 

4.1 A brief presentation of existing formalisms for the description of particle oscil- 
lations 

4.1.1 The plane v^ave description of particle oscillations 

We present the elementary treatment of particle oscillations based on plane waves, which is the most well- 
known and widely used due to its simplicity. This approach is known to be fundamentally inconsistent [6l 

13 



160] . as it relies on infinitely extended wave packets with zero momentum spread in each mass eigenspace. 
The wave packet methods, presented in Sec. 4.1.2, and the quantum field theory methods, presented in 
Sec. 4.1.3, were developed partly in order to resolve outstanding issues in the plane wave treatment and 
to place the derivation of the oscillation formulas on a firmer mathematical ground. 

By definition, an oscillating particle is described in terms of a Hilbert space % = (BiT-ii, where the 
subspaces Tii carry irreducible unitary representation of the Poincare group with the same spin but 
different masses mi. Ignoring for simplicity the spin degrees of freedom and particle decay, each mass 



eigenspace Tii is spanned by the generalized eigenstates \p,i) of the Hamiltonian Hi = JfP' + m? for a 
free particle of mass rrij. One defines 'flavor eigenstates' 

|p,a)=5]C/*,|p,i), (18) 

i 

where Uai is the mixing matrix appearing in the interaction Hamiltonian for the oscillating particles. 
The transition amplitude between a generalized eigenstate of flavor a ,at time t = and location x = 0, 
and a generalized eigenstate of flavor /3, at time t and location x = L, is 

^,^^(t,L) = (p,/3|e-*^*+^P-i^|p,a) = e^P'^ ^ e-^^'*C/^C/^,(p,i|p,i), (19) 



where Ei = Wp^ + m?. The term (p, i|p, i) is really infinite, but in the plane wave method we are 
interested primarily in the relative phases between the mass eigenstates, so this issue can be ignored. 
The squared amplitude 

|^„^^(t,L)|2 = 5]C/*,C/^i[/,,[/^^.(p,i|p,i)(p,j|p, j)e-^^'*+^^^* (20) 

depends explicitly only on t. It would also carry a dependence on L if we had assumed different values of 
Pj in the definition of fiavor eigenstates, Eq. (jlSp . but this is not necessary for the present demonstration. 
In order to relate the square amplitude Eq. (I20p with a measurable quantity, it is necessary to convert 
its time dependence into spatial dependence. There is no rule of quantum mechanics that allows one to 
do so, so one has to rely on ad hoc prescriptions. To this end, we note that each subspace is characterized 
by a velocity Vi = p/Ei, where p = |p|. At the classical level, this corresponds to the time of arrival 
ti = L/vi, for a particle starting at x = to a detector at x = L, where L = |L|. 

The standard oscillation phase is obtained by substituting the time parameter t with a mean time 
of arrival i at the detector, e.g., for the time of arrival calculated for a mean energy E, — this procedure 
is usually called an 'equal tiines prescription'. We then obtain a phase difference 

^i - ^j = [Ei - Ej)i = {Ei - Ej)LE/p, (21) 

which renders the square amplitude Eq. ()20p an oscillating function of L, from which the standard 
oscillation wave-number kij = —\Ei — Ej\ follows. In the ultra-relativistic regime, relevant to neutrinos, 
the oscillation wave-number is 

\mf-mj\ 
kij = — . (22) 

On the other hand, if we substitute the different values ti = L/Ei for arrival time in each mass 
eigenspace (unequal-times prescription), one obtains (j)i — (j)j = -Ejtj — Ejtj = {Ef — E'^)L/p = {mf — 
m'^)L/p, corresponding to an oscillation wavenumber 



I 2 2i 
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which is twice the standard expression Eq. ()22p in the ultra-relativistic hmit. 

We note that the non-standard result Eq. (j23p includes the contribution from interferences in the 
time of arrival corresponding to the different values ti in the different mass eigenstates. The derivation 
of the standard result Eq. ([22]) involves an assumption that such interferences are somehow suppressed. 

4.1.2 Wave-packet methods 

The natural generalization of the plane wave method is the consideration of amplitudes between properly 
square-integrable wave-packets (BJ [Ml EH ESI EHj . That is, one defines Hilbert space vectors 

|*^(0),a) =5:C/^*f(x)|x,i), |M/^(L)/3) = 5:C/^*f(x)|x,i), (24) 

i i 

in terms of profile functions "^^ and ^^ centered around the production ( x = 0) and detection (x = L) 
region respectively; \x,i) stand for the generalized eigenstates of position in the i-th mass eigenspace. 
One then computes the amplitude 

Aa^ls{t,L) = (M'^(L),/3|e-^^*|M'^(0),a). (25) 

The simplest case that allows for analytic calculations corresponds to isotropic Gaussian wave-packets 

vl'f (x) = -^e ^ , *f (x) = -^e ^^ , (26) 

where one assumes that the production region is centered around x = with spread equal to ap and 
the detection region is centered around x = L with spread equal to an- 

The amplitude Eq. (I25p can be computed in the no-dispersion regime, defined by the approximation 
Ei{p) ~ Ei{pi) + Vi ■ (p — Pi). Restricting for simplicity to one spatial dimension, one finds 



A^Mt,L) = y^5:^,,f/^e-^»*^---^^^, (27) 

i 

where a"^ = aj^ + dp. 

In order to express the probabilities in terms of observable quantities alone, one needs to remove the 
time-dependence from the amplitude in Eq. (I27p . This can be done in two ways. 

1. The interference phases in the modulus square of the amplitude Eq. (j27p 

^i - ^j = {E, - Ej)t - L{p, - pj). (28) 

are the same with the ones obtained from the plane wave methods. They can, therefore, be 
evaluated with the same methods, i.e., by making a specific choice of time t corresponding to a 
'mean' time of arrival. The standard expression for the oscillation wavelength necessitates the 
postulate of an additional prescription for the terms in Eq. (j28p . The most usual form of such 
prescription is the assumption of equal-times ti = tj, i.e. an assumption that interference takes 
place only for equal propagation times J 12l . Alternatively, the standard oscillation formula may be 
obtained by an equal energy assumption Ei = Ej [IHIES], by an equal velocity assumption Vi = Vj 
|67j . or by a specific choice of 'mean velocity' for the wave packet [151 [161 184] . 

2. The introduction of ad hoc prescriptions is avoided in the wave-packet method by integrating the 
squared amplitude \Aa^/3{t, L)\'^ over time fBB], i.e., by defining a probability for particle detection 

Pa^piL) = C f dt\Aa^pit,L)\\ (29) 

Jo 

in terms of some long integration time T. The constant C is introduced for purposes of normal- 
ization of the probability to unity. 
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4.1.3 Quantum field theoretic methods 

The oscillating particles are not directly observable, therefore, the natural framework for the treatment 
of particle oscihations is quantum field theory [MlEOlEIlEZlESlIHESESlEaEHlEalEOllSOlllS]. In 
the latter, oscillating particles are described to internal lines of Feynman diagrams, whose external lines 
correspond to incoming particles in the production region and outcoming particles at the detector. 
In this approach, one computes amplitudes of the form 

{^°,tD\Te-'I'^'''^'\^^,tp), (30) 

where |^^,tp) describes the state of the incoming particles in the production region, |^^,i£)) describes 
the state of the particles produced in the detection process, T denotes time-ordering and Tij is the 
interaction Hamiltonian for the total process. Typically, |^^,tp) correspond to wave-packets localized 
at the production region at a specific time tp and I^Di^d) corresponds to wave-packets localized at the 
detector region at a time tp). The calculation of the amplitude Eq. ([30]) to leading order in perturbation 
theory introduces the characteristic dependence of the amplitude on the mixing matrix Uai in Eq. (I27p . 
Usually, one constructs the relevant probabilities from an integration of the amplitudes Eq. ([30]) squared 
over time T = tp, — tp. 

Ref. [To] avoids the use of integration over time in squared amplitudes by following a different 
procedure. One computes the S- matrix elements for the weak process n -^ p + e~ + !>, for the production 
of the antineutrino localized around x = and the subsequent scattering process v' + e~ -^ u + e~ 
localized around x = L for the antineutrino detection without reference to an initial and final time. The 
S'-matrix formalism in this case implies that the initial and final states are 'in' and 'out' respectively, 
i.e., they are defined at i — )• — oo and t — )• +oo. The authors then construct the scattering cross-section 
at the limit L — ?■ oo recovering the standard result. Ref. [75j avoids the S'-matrix formalism, relying 
instead on time-dependent perturbation theory. 

The quantum field theoretic method provides a substantial improvement in the theoretical treatment 
of particle oscillation, because one needs not introduce ill-defined flavor eigenstates and also because 
it allows for a more elaborate treatment of the production and detection of oscillating particles — for 
a comparison of the methods, see Ref. [62]. A quantum field theoretic treatment is essential to a 
first-principles treatment of the problem and we shall employ it in Sec. 5.2. 

However, there is little difference at the abstract level of formalism between the quantum field the- 
oretic and the wave packet methods. Both make the same use of vector states and evolution operators 
in the definition of probabilities, essentially presupposing a description in terms of von Neumann mea- 
surements. This elementary point may be obscured by the S matrix language commonly employed in 
quantum field theory, so we elaborate here on the meaning of the amplitude Eq. (]30p . In quantum 
theory, the dependence of a vector state on time arises from unitary time evolution. The vector states 
|^D,P)i_D,p) are then states of the form e~*''*|^^p) for some vector states |^^p) and a Hamiltonian 
h. The choices for \^D,p,tD,p) in the literature imply an identification of h with the Hamiltonian Hq of 
free particles. Then |^^p) corresponds to single-particle states for outcoming and incoming particles, 
centered in the detection and production region respectively. The amplitude Eq. (j30p is then equivalent 
to 

(^'o|e-^^«*°re-*/'^*^^e^^«*^|^'p), (31) 

where Hj = J d^xTii. 

Eq. ()3ip is a meaningful quantum expression only if the integration in the time-ordered exponential 
is for the time interval [tp,tD], in which case, Eq. ([3T]) equals (^'p,|e~*(-'^°"'"-'^^)(*^~*^)|^p). The squared 
amplitude then provides the usual expression for probabilities in von Neumann measurements. Hence, 
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as far as the description of measurements is concerned, both wave packet and quantum field theoretic 
methods employ the same assumptions. 

4.2 Critical analysis of the existing formalisms for particle oscillations 

4.2.1 Justification in terms of first principles. 

The first point of critique towards the particle-oscillations methods described above is that the construc- 
tion of probabilities does not follow from the first principles of quantum theory. The vast majority of 
these methods identifies the relevant probabilities with the modulus square of amplitudes, such as Eq. 
()25p . However, this expression is relevant only for measurements that take place at a pre-specified mo- 
ment of time (or time-interval) — see. Sec. 3.1. In particle oscillations, the experimentally fixed quantity 
is the location of the detector, not the time that a particle is detected. The time of detection is a variable 
rather than an external parameter. This inconsistency in the definition of probabilities is the reason why, 
in order to obtain a connection with measurable quantities, one needs to remove the time dependence 
from the square amplitudes — either by introducing additional prescriptions or by time integration. 

The construction of probabilities using squared amplitudes is the most common practice in applica- 
tions of quantum theory, but it is neither necessary nor unique. The most general probability assignments 
in quantum theory correspond to POVMs (see. Sec. 3.1.3), and these can be applied to cases where 
the simple definition in terms of squared amplitudes fails. We maintain that the correct probability 
assignment for each particular case should be constructed from first principles, in order to reflect pre- 
cisely the physics of the experiment. The use of a probability assignment without prior justiflcation may 
lead to inconsistencies in the correspondence of elements of the formalism to experimentally determined 
variables. Hence, the modeling of particle oscillation experiments by von Neumann measurements is 
problematic, because the setup of the former violates a main assumption in the latter, namely, that 
measurements take place at a pre-determined instant of time. 

Even should one attempt to model particle oscillations by von Neumann measurements, ampli- 
tudes such as Eq. ([251) or Eq. (pOj) are insufficient: they ignore the effect of the detector, which 
is a macroscopic system with a large number of degrees of freedom. To see this, we note that the 
squared amplitude {Aa^pit, L)\'^ can be brought in the standard form Tr{pIlL) by identifying the state 
p = e~*^*|*'^(0))(^^(0)|e*^* and the projector Ul = \'^^{L)){^^{L)\, that describes the measurement 
scheme. The projector H/, is one-dimensional {Trtli = 1), i-e., it corresponds to an extremely fine 
measurement that determines a phase space volume for the particle of size h [18^ [29] . One-dimensional 
projectors may be relevant for sharp measurements (observables with discrete spectrum), but they are 
unsuitable for measurements of continuous variables. As explained in Sec. 3.1, the interaction with 
the measuring apparatus necessitates the use of much coarser projectors (TrH >> 1) for continuous 
variables. Hence, the use of the squared amplitude |^a_>^(t,-L)p for the probabilities misrepresents the 
coarse-graining inherent in the detection process. 

For example, a projector describing an unsharp measurement of position is 

Hl = / dyfiy - L)|M/^(y))(vI/^(y)|, (32) 

where /(y) is a positive function centered around y = with a spread of order 6L [81j. The corresponding 
probabilities are 

TriptlL) = J dyfiy - L)\Aa-,^{t,y)\^. (33) 

Such modifications in the derivation of probabilities would not affect the phases relevant to the oscillations 
(as long as 5L « L), but it would make a difference in the expression for the so called 'coherence length' 
[6l [631 ESl ESI [82] i-e-, a scale L^oh such that, if kijLcoh >> 1 no oscillations are observable. 
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4.2.2 Times of arrival and the different prescriptions 

In Sec. 4.2.1, we argued that the definition of probabihties from the squared ampHtude \Aa^(3{t, L)\'^ 
is not suitable for particle oscillation experiments. Its use presents an immediate problem of reconciling 
the theoretical description with the experimental procedure. In particular, it is necessary to remove 
the time dependence from the squared amplitudes \Aa^i3{t, L)\'^ ■ In this section, we examine the most 
common procedure employed for this purpose, namely, the substitution of the time parameter t in the 
squared amplitude by the 'mean' or classical time of arrival for a particle propagating from the source 
to the detector. The second common procedure of eliminating the time variable, by time-integration, is 
considered in Sec. 4.2.3. 

The substitution of a mean time of arrival in the squared amplitudes \Aa^i3{t, L)\'^ has no funda- 
mental justification. In fact, this procedure is inconsistent with quantum theory. The time of arrival 
depends on the initial state (for example, its mean energy), and its substitution into the squared am- 
plitudes lAa^pit, L)\'^ leads to a probability assignment that is not a linear functional of the system's 
density matrix, thus contradicting a fundamental requirement for any probability assignment in quan- 
tum theory. In particular, as we explain in more detail in Sec. 4.2.3, a probability assignment that does 
not respect linearity does not have a consistent statistical interpretation. 

The lack of a solid mathematical justification for the aforementioned procedure implies that there is 
no independent theoretical criterion to settle any ambiguities that may arise in its application. This is a 
problem, because, as we saw in Sec. 4.1.1, a variation of this procedure that is equally plausible a priori — 
the so-called unequal-times prescription — leads to a prediction for the particle oscillation wavelength that 
differs from the standard one by a factor of two. Hence, a significant part of the related bibliography 
is devoted in showing that the unequal-times prescription cannot be valid, aiming to strengthen the 
standard result by the exclusion of its alternative. 

In what follows, we will show that (i) that the physical arguments developed against the non-standard 
oscillation formulas do not apply to our result which is obtained by a different method, (ii) the standard 
derivations ignore the probabilistic and quantum aspects of the time of arrival, and (iii) the timescales 
involved in particle oscillation experiments suggest that interferences in the time of arrival contribute to 
the oscillation wavelength thus leading to a non-standard particle oscillation formula. 

A strong argument against the unequal-times prescription is that it leads to an unphysical prediction 
of oscillations for recoil particles; for instance, the muon in the reaction vr — t- /ii^ or the A in the reaction 
ir~p — )• Ai^". This is highly counterintuitive, and it is generally considered as a severe deficiency of the 
method. The formalism we develop here leads to a non-standard oscillation formula, but it does not 
predict oscillations for recoil particles — see Sec. 5.3. A careful description of the measurement process, in 
terms of particles that are actually being detected, shows that oscillations are only possible for particles 
with different mass eigenstates. Therefore, we maintain that the oscillations for recoil particles is the 
artefact of an improper mathematical framework for the evaluation of quantum probabilities. It is not 
a consequence of the general statement that different values of the time of arrival contribute coherently 
into the particle oscillation probability. 

Common arguments against the unequal prescription involve statements that quantum interference 
only occurs between states defined at the same spacetime point [121 113^ 183]. or that measurements take 
place in space and not in time [66j. However, the restrictions implied by such statements do not follow 
from the basic principles of quantum theory. Any observable is recorded by pointer variables on a 
measuring apparatus, and interferences are observed through the distribution of values of these variables 
in a large number of runs. Nothing in quantum theory forbids a pointer variable from becoming correlated 
with properties of the measured system that correspond to different instants of time. As shown in Sec. 
3.2.2, the time of arrival may exhibit interferences, just as any other quantum observable. In general. 
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we maintain that any statements about what cannot be measured in quantum theory can be made only 
on the basis of a precise modehng of the measurement process. 

In Sec. 2.1, we stated that a key problem in existing approaches to particle oscillations is that they 
conflate the time-parameter of Schrodinger's equation with the time of arrival, which is a physical observ- 
able. In the present context, this problem is manifested by the ad hoc substitution of the time parameter 
t — arising from the solution of Schrodinger's equation — in the squared amplitudes {Aa^pit, L)\'^ by a 
'mean' time of arrival. This procedure ignores the fact that the time of arrival is a genuine quantum 
observable. Indeed, the time of arrival is determined by the coincidence of a microscopic quantum event 
(for example, a particle reaction in the detector) with the reading of an external clock. While the value 
of the time of arrival is a classical variable, it is correlated to a microscopic process that is governed by 
the rules of quantum theory. Consequently, the values of the time of arrival exhibit fluctuations, and 
these must be described by quantum probabilities. The substitution of a mean value for the time of 
arrival in the squared amplitudes \Aa^(3{t, L)p ignores all effects pertaining to quantum fluctuations; in 
particular, it excludes a priori the possibility of quantum interference in the values of the time of arrival. 

Acknowledging the time of arrival as a quantum observable, brings new light in the distinction 
between the standard oscillation formula and the non-standard one following from the unequal-times 
prescription. The former essentially assumes that quantum interferences in the time of arrival are either 
suppressed or negligible; the latter implies that such interferences persist and they contribute significantly 
to the particle oscillation probability. As with any other quantum variable, the issue whether such 
interferences persist or not can be resolved only through a careful study of the physical situation at 
hand, namely, of the quantum system and the measuring apparatus. To this end, a fully quantum 
description of the time of arrival is necessary. This, we undertake in Sec. 5. 

An additional problem in the use of a 'mean' time of arrival in the squared amplitudes \Aa^(3{t, L)\'^ 
is that the very notion of a 'mean' time of arrival is ambiguous in absence of a full quantum treatment. 
For example, in Refs. [151 [161 [H] the standard expression for particle oscillation wavelength is obtained 
by substituting a mean time of arrival i = L/v in Eq. (|20]) . where L is the source-detector distance, and 
the mean velocity v of the wave-packet is defined as 

El + E2 

where the indices 1 and 2 refer to two mass eigenspaces. Eq. (|34p has no fundamental justification 
in quantum theory, where the only notion of a mean velocity is the expectation value of the operator 
pH~^. There is no reason why the mean velocity should not be defined by any other expression, such 
as the the arithmetic or the geometric mean of the velocities Vi = Pi/Ei. A justification for Eq. (j34p is 
that it guarantees that no oscillations for recoil particles are predicted, without making an equal-times 
assumption. If, however, oscillations of recoil particles is an artefact of a mathematically ambiguous 
formalism, as our results in Sec. 5 indicate, this argument loses much of its power. 

In our opinion, one of the most important arguments in support of prescriptions for the standard 
oscillation formula (it is also quoted as definitive in the Review of Particle Physics [5j) is essentially 
the following [T6]. Consider Eq. (j28p for the oscillation phases for two mass eigenspaces, i,j = 1,2. 
Substitute for i the expression 

i=- + 5t, (35) 

V 

where v is given by Eq. (|34p and 5t « L/v is "a correction to describe any discrepancies, including 
quantum fluctuations" . Then, 

rr^-ml Ef - EJ 

"^1 " ^^2 = ———L - -^r-—ErSt- (36) 

Pi +P2 -C/1 + -C/2 
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The first term in the right-hand-side of Eq. (j36p is the standard oscillation phase, while the second term 
is a correction much smaller than the first term. The correction 6t is deemed to be also much larger than 
the time-scale 6ti2 = \L/vi — L/v2\, that corresponds to the different centers of the wave-packets in the 
different mass eigenspaces. The term 6ti2 is responsible for the different prediction of the unequal-times 
prescription. One then concludes that the effect of 6ti2 in the oscillation phase is negligible and thus 
the standard oscillation formula holds. 

We note that the introduction of the correction 5t in the phase difference is essential for the argument 
to work. However, if St corresponds to the fluctuations of the time of arrival, its introduction into the 
phase difference is not justified. For the argument's sake, let us assume that the parameter t in the 
amplitude Aa^i3{t, L) can indeed be substituted by the time of arrival. Let p{t) be the probability 
distribution describing any fluctuations of the time of arrival — 6t corresponds to the mean deviation 
of this distribution. According to standard probability theory, the incorporation of these fluctuations 
involves weighting the squared amplitudes \Aa^(s{t, L)\'^ with p{t), i.e., one should consider the integral 

dtp{t)\Aa-,p{t,L)\\ (37) 

Then, the integration Eq. (I37p does not lead to the introduction of the mean deviation 5t into the phase 
difference (pi — (p2- For example, if p{t) is a Gaussian ~ exp[—{t — L/v)'^ /(26t'^)], the integration in Eq. 
()37p would only lead to a suppression of the oscillation's amplitude. Hence, the fact that St « L/v or 
that 5t >> 5ti2 has nothing to do with the relative size of terms in the oscillation phase, and, hence, 
with the oscillation wavelength. 

Moreover, the specification of a single parameter 6t for the description of the fluctuations in the 
time of arrival is insufficient if these fluctuations are quantum. A single parameter does not allow 
for distinction between incoherent and coherent fluctuations because these have very different physical 
behavioo The former tend to suppress oscillations while the latter would contribute additional terms to 
the oscillation phase. Ref. [16] emphasizes that the understanding of time fluctuations is fundamental for 
the proper theoretical description of particle oscillations. We strongly agree on this point. Furthermore, 
we maintain that such a treatment should follow precise probabilistic reasoning and it should take into 
account the quantum nature of the time of arrival. 

The treatment of the time of arrival as a physical observable has important implications for the 
description of particle oscillations. A particle-detection event takes place at a specific locus and a 
specific moment of time. In the context of particle oscillations the detection time corresponds to the 
time of arrival. Since the latter is a quantum observable, the detection process corresponds to a joint 
measurement of time of arrival and position, which can be described in terms of a suitable POVM — 
see. Sec. 3.1.3. For a joint measurement, the separation of the timescales 6t and 6ti2, noted in Ref. 
|16j . leads to the opposite conclusion from the one reached in that reference. If the sampling of one 
of the two variables in a joint measurement is much coarser than the relevant interference scales in 
the quantum state, then the measurement of this variable does not affect the probability distribution 
for the other variable. In other words, the fact that the timescale St and 6ti2, defined earlier, satisfy 
dt » 6ti2, implies that the measurement of the time of arrival does not disturb the contribution of the 
time-of-arrival interferences in the oscillation probability. 

To see this, one may consider the two-slit experiment. Let 5L be the distance between the two slits. 
Consider now a special case of a joint measurement: first, a recording of the position of a particle as it 
passes through the slits and its later position on a screen. If the accuracy d of the first measurement 
satisfies d » 6L, the slit through which the particle crossed is undetermined, so the interference 



^For example, a Gaussian wave function exp{—x^/L^) has approximately the same spread L with the superposition 
exp{—x^/a^) + exp{—{x — L)^ /a^), if a « L, but their physical properties are very different. 
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pattern in the second measurement will persist. Only for a finer-grained first measurement d ^ 5L are 
the interferences suppresseco- The key point in this version of the two-slit experiment, is that we also 
have a joint measurement of two observables. A coarse measurement of one observable does not affect 
the coherence in the measurement of the second observable. 

4.2.3 Integration of probabilities over time 

In Sec. 4.2.2, we examined the procedure of removing the time-dependence from the squared amplitude 
\Aa^l3{t, L)\'^ by substitution with a 'mean' value for the time of arrival. In this section we consider the 
alternative procedure of integrating the squared amplitude \Aa^i3{t, L)p over time, as in Eq. (f29]l . This 
procedure is more common in the recent works on particle oscillations, and it is often applied in a field 
theory setting. The arguments presented here do not depend on the way such amplitudes are calculated, 
but follow from basic properties of the formalism of quantum theory. 

From a mathematical point of view, \Aa^i3{t,L)\'^ is not a density with respect to t, so integration 
with respect to t is ill-defined. This problem arises because t in \Aa^i3{t, L)\'^ is an external parameter and 
not a random variable (i.e, an observable). The integration over parameters (rather than over random 
variables) does not define proper probability densities, neither in classical nor in quantum probability 
theory. The reason for this is that the probabilities for different alternatives can be added only if the 
alternatives are disjoint. If the alternatives are not disjoint, the joint probability is not obtained by the 
sum of the individual probabilities 

An elementary example from ordinary probability theory serves to illustrate that events distinguished 
only by the value of the time parameter t fail to be independent. Let us consider a variable X taking 
values in some set F and consider probabilities defined on paths X{-) of such variables labeled by the 
time parameter t. The event (alternative) "X = x at time ti" is not disjoint to the event that "X = x at 
time t2" • There are paths that satisfy both propositions above, for example a path that satisfies X = x 
at all times t. One cannot therefore obtain the joint probability for the event "X = x at ti or X = x at 
time ^2" by summing the individual probabilities. 

The same also holds in quantum theory [86j. As a result, the integral in Eq. ([29]) does not lead 
to well-defined probabilities. In particular, it does not preserve normalization: to normalize one has to 
divide by a constant C. In general, the normalization constant C depends on the initial state of the 
system, so the probability assignment is a non-linear functional of the initial density matrix p, a property 
that is incompatible with the rules of standard quantum theorjo. 

Linearity with respect to the initial state is an essential property for any probability assignment. It 
guarantees that one can combine sub-ensembles through statistical mixing. This means the following: let 
p^ be the density matrix corresponding to a sub-ensemble of events distinguished by the macroscopically 
determined quantities ^, and let Pp^{A) be the probability distribution corresponding to an observable 
A — see the Appendix B for details. The ensemble constructed by weighting each sub-ensemble with a 
factor w^ is described by the convex combination J2^w^p^. Linearity guarantees that the probability 



^This property is seen most clearly in the so-called "which- way" experiments [SSJ, where a microwave cavity is used in 
order to perform a measurement of position as the particle crosses through the slits: the state of the electromagnetic field 
in the cavity determines the accuracy of this measurement. 

■^For two overlapping alternatives Ui and U2 and probability distribution p(), the joint probability p{Ui U U2) = 

p{Ui)+p{U2)-p{UinU2). 

^One might respond to our argument by stating that the normalized probabilities Eq. (|29|) are in fact conditional 
probabilities at time t, provided that detection has taken place. But in this case one should prove that the normalization 
constant C is proportional to the total detection probability, as predicted by quantum theory. This is problematic because 
C has no analogue in standard expressions for detection probabilities. Moreover, the recourse to conditional probabilities 
does not remove the problem of defining unconditioned probabilities. In particular, a proper normalization of uncondi- 
tioned probabilities is necessary for the description of experiments at which neutrino detection takes place at two different 
locations — see, for example, the proposal in Ref. [87) . 

21 



for the variable A in the total ensemble will also be a convex combination ^ t "^^Ppe (^) • The lack of 
linearity in the probability assignment Eq. (j29p implies the lack of a consistency in the combination of 
different subsets of data [88j . 

The requirement of linearity and positivity in the probability assignment is very restrictive. It implies 
that any probability assignment must correspond to a POVM. In the present context, the POVM should 
describe a joint measurement of position L and arrival time t. Given such a POVM n(t, L) we can 
integrate consistently over time of arrival in the probabilities Tr [/5II(t, L)] obtaining a probability density 
that depends only on position 

p{L) = f dt Tr[pl[{t, L)]. (38) 

Jo 

It is only in the context of equations of the form Eq. (|38p that integration of probabilities over time can 
be justified in accordance with the rules of quantum theory. 

The problems above would be resolved if one constructed amplitudes that depend on the time of 
arrival rather than the parameter t of Schrodinger equation. However, neither in this case would Eq. 
()29p be justified. In quantum theory a summation of different alternatives at the level of probabilities is 
justified only when the alternatives are macroscopically distinct (for example, if they are distinguished 
by a measurement scheme). In absence of such distinction, the summation ought to take place at the 
level of amplitudes. This property is clearly manifested in the two-slit experiment. If we can determine 
macroscopically through which slit the particle passed, we sum over the corresponding probabilities. 
If we cannot, then we sum over amplitudes. In the latter case, we observe an interference pattern, in 
the former we do not. Hence, the assumption of an integral over time, at the level of probabilities, 
presupposes that we can distinguish temporal alternatives at a scale smaller than any interference terms 
in the system's state. This is not the case in particle oscillation experiments. As also discussed in 
Sec. 4.2.2, the temporal resolution of any particle detection corresponds to a time-scale r^es which is at 
best of the order of lO^^s for a solid state detector. On the other hand, the characteristic timescale of 
interferences 6tij is of the order of W~^^s, i.e., many order of magnitudes smaller than r. 

Since the temporal resolution does not suffice to decohere the interferences, quantum theory would 
suggest an integration over different temporal alternatives at the level of amplitudes. Indeed, in this 
paper we describe the probabilities relevant to particle oscillations through a POVM II(t, L), and we 
show (see Sec. 5. 2) that in the physically relevant regime, the resulting probabilities are proportional to 
I /q dtAa^idit, L)\'^. a non-standard expression for the particle oscillation wavelength then follows. 

Finally, we comment on a recent work Ref. ^62j, which contains a very thorough treatment of the 
normalization in the time-integrated probabilities. In that work, a normalized probability distribution 
is obtained by integrating squared S-matrix amplitudes with respect to both production time tp and 
detection time to- Here, one assumes that the particle flux is constant at both the production and 
the detection region. This methodology applies in the ultra-relativistic regime and for quasi-degenerate 
neutrino masses; in the general case, the oscillation probability is undefined. 

In our opinion, the assumption of a constant particle flux at detection is too strong. Such an 
assumption is reasonable at a highly coarse-grained, macroscopic time-scale but there is no justification 
in its extrapolation to small time scales. In particular, this assumption excludes a priori possible 
variations of the flux resulting from different 'propagation velocities' in different mass eigenspaces, i.e., 
the potential contribution of interferences in the time of arrival. 

4.2.4 The inapplicability of the S-matrix formalism 

In quantum field theoretic treatments of particle oscillation, one often applies the S'-matrix formalism. 
The S-matrix is well-suited for the construction of probabilities in scattering experiments. Its domain 
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of validity, however, is not universal: it has been constructed for the description of ordinary scattering 
experiments in which the detector is far from the scattering region, the measured observables are the 
outcoming-particles' momenta, and one is not interested in the location and timing of the scattering 
event. In contrast, in particle oscillation experiments the main observable is the location and timing of 
the reaction through which the particle is detected. 

To this end, we consider a simple model for the detection of particles produced from the scattering 
of a beam into a target, using the formalism of Sec. 3.1. In this case, the particle detector is located at 
a macroscopic distance from the scattering region. Let to be an arbitrary time after the scattering has 
taken place, so that after tg evolution is governed by the Hamiltonian Hq of free particles. 

We assume that the measured observable A is momentum (which commutes with the free particle 
Hamiltonian), or indeed any other quantity such that [A, i^o] = 0. Since the time-of-arrival is not a 
measurable quantity, the description of measurement in term of a profile function f{t), as in Sec. 3.1, 
suffices. We assume that the function /(t) takes non-zero values only for t > Iq, and that t = T is 
an upper limit to the time the particles spend in the detection region. It is easy to show that, since 
[A, H] = 0, Eq. ([II]) becomes 

p{X) = TrHs {puMX - bA)) , (39) 

i.e., the values of A are correlated with the values of the rescaled pointer observable X/b, where b = 
J dtf{t). The density matrix pt^ in Eq. (j39p arises from the evolution of the initially prepared state 
under the full Hamiltonian H including interaction terms: pt^ = e~*^*°/5oe*^*''. The key point in Eq. 
()39p is that there is no dependence of the probabilities on an instant of measurement or on the duration 
of the measurement. Since the observable A commutes with the Hamiltonian, the measurement records 
the distribution of A^ contained in the state pig, after scattering. Since to is itself arbitrary, we can make 
it arbitrarily large, so that pt = SpqS\ where S is the 5-matrix. 

This argument perhaps belabors the obvious, but it demonstrates that 5-matrix formalism provides 
accurate probabilities without any reference to an instant of measurement, only for observables that 
commute with the free-particle Hamiltonian. For any other observables (in particular, position), the 
probabilities remain dependent on the time the measurement takes place and the limit t — >• oo is unwar- 
ranted. Therefore, the application of the S matrix formalism in the context of particle oscillations, where 
the relevant observable is the location of the scattering event, takes it outside its domain of validity. For 
another critique of the use of the S-matrix, see Ref . [89J . 

Another way to see the unsuitability of the 5-matrix for the description of particle oscillations follows 
from Eq. ()3ip . This equation is quantum mechanically meaningful only if the domain of integration in 
the exponential is \tp,tD\- However, the time integral in the field theory treatments extends along the 
full real axis, i.e., the S matrix formalism is employed. In the perturbation theory treatment of scattering 
experiments, the integral j]-^ dt in the time-ordered exponential is substituted for /^ dt because one 
assumes that ti and tj are much earlier and much later respectively than the time of the scattering event, 

-i r*^ Hjdt 

SO that the interaction terms are effectively "switched-off'. In this case, substituting Te *» with 

the S'-matrix is an excellent approximation. However, in a particle oscillation experiment tp and tp, are 
the times at which scattering events are assumed to take place. There is, therefore, no justification in 

^ -« r*° ■f^/'^* ■ , , rr 

exchanging Te *p with the o-matrix. 

The issue of the applicability of the S matrix formalism to particle oscillation experiments is a 
key problem in quantum field theory treatments that do not employ integration over time in their 
derivation [70^ I73j . In particular, the interpretation of the 5-matrix elements on spatially localized 
states as corresponding to position measurements takes the 5-matrix formalism outside its domain of 
applicability. 
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4.2.5 Special initial states 

A class of arguments in favor of an equal-energy prescription |131 190] states that the time of emission of 
a neutrino is unknown, so one should average the initial state over time. This leads to an initial state 
that is a sum over incoherent energy eigenstates, so that the energies of the initial state are the same in 
all mass eigenspaces: Ei = Ej. We find the following problems with this argument. 

First, in any particle production process we always obtain some information, such as macroscopic 
records of position and momentum of other particles produced in the reaction. In the Appendix B, 
where we discuss in detail the initial state issue, we denote these parameters collectively as ^. Then, 
for each value of ^ there is a different density matrix p^ corresponding to the sub-ensemble of neutrino 
creation events characterized by macroscopic records ^. The total ensemble is then a convex combination 
Sf P^P^i 1 > Pg > 0. There is no general argument why the expression J^^P^P^i referring in general to a 
large number of observed values ^, and corresponding to different particle production processes, would 
always be proportional to an integral / dtp{t) over the initial state of the particles. 

Moreover, averaging a density matrix over the time of particle production suffers from the same 
mathematical and conceptual problems as averaging the density matrix over time of arrival; the dis- 
cussion in Sec. 4.2.3 applies to this case as well. Like the time of arrival, the time of production is a 
physical observable and it does not coincide with the time parameter of Schrodinger's equation. A proper 
description of the dependence of probabilities on the time of particle production to should proceed from 
the inclusion of to into the variables ^ of the POVM H^, that corresponds to the variables determined 
during the production process. Only with such a treatment would it make sense to average the state over 
the production time Iq, but in this case, integration over io would be a part of the convex combination 
defining the neutrino initial state. 

In general, the recourse to a fine specification of the initial state is undesirable in any theoretical 
prediction or explanation, because it involves strong assumptions that cannot be independently verified. 
It is highly unlikely that such fine-tuning could persist in different experiments characterized by different 
preparations of the system. Even the justification of the properties of the neutrino wave-packet by the 
physics of the production process [78], requires the introduction of a number of unobserable parameters 
and a choice of specific regimes. In the method we develop in this paper, the results do not require 
any fine-tuning in the parameters of the initial state. We obtain the non-standard oscillation formula 
even if we assume an initial state for the neutrinos with equal energies Ei in different mass-eigenspaces. 
In fact, our results are valid for a large class of initial states that only share some broad qualitative 
characteristics. 

4.2.6 Time integration in probability currents. 

A different class of theoretical approaches for the description of particle oscillations considers the particle 
flux at the detector, defined by time-integration of a probability current J^(,x, t) for flavor a — see, Ref. 
|91j . The current JJy{,x,t) for flavor a is constructed as to (approximately) satisfy a continuity equation 

^^^ + V.j;(x,t) = 0, (40) 

where pa{^,t) is the quantum probability distribution at time t. J^ is interpreted as a probability 
current. One then defines a flux for the particles of flavor a, detected in a region A, at distance L from 
a source, as 

^a{dA)= [ dt f d\^Jl{x,t), (41) 

JO JdA 

where T is an integration time for the experiment and dA the boundary of the detector region A. 
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There are open theoretical issues in this approach: (i) how to properly define flavor currents in the 
relativistic regime, relevant to neutrinos [92], and (ii) how to interpret the fact that the conservation 
equation is not exact. The latter issue, in particular, creates problems in the normalization, since the 
total flux along a closed surface S around the production region is not conserved. 

However, even if such problems were resolved, a more fundamental problem in the consistency of the 
integration over time in Eq. ()4ip would have to be addressed. Probability currents in quantum theory 
may also take negative values [93], so the flux in Eq. ([^T|) does not, in general, correspond to the total 
number of detected particles. This issue has been discussed in relation to the definition of probabilities for 
the time of arrival |94il951 [96]. In that context, it has been argued that in regimes at which probability 
currents are approximately positive, they do yield physically acceptable distributions for the time of 
arrival ^97j. This procedure, however, cannot provide a rule for computing the expected particle number 
for a general initial state. Moreover, the probability current defines positive probabilities only in the 
quasi-classical regime for the time of arrival [55]. In contrast, particle oscillations are a consequence of 
quantum interference, and one should expect that they cannot be adequately described in the quasi- 
classical regime. 

4.2.7 Quantum models of the production and detection process 

A variation of the quantum field theory approach to particle oscillation involves the description of the 
production and detection events in terms of quantum systems playing the role of 'sources' and 'sinks'. 
These can be modeled, for example, by a pair of families of harmonic oscillators, labeled by flavor, and 
localized in the production and detection region. The interactions of the oscillators with the neutrino 
field is switched on by profile functions dependent on time |72j . The coupling of the oscillators to the 
neutrino field is constructed so that the desired dependence of the probabilities on the mixing matrices 
Uai is obtained. We denote by |na,?n,^,0) the state for the source oscillator for the a flavor in the n-th 
energy eigenstate, for the sink oscillator for flavor /3 in the m-th energy eigenstate and for the neutrino 
vacuum, one constructs the amplitudes 

^„^^ = (l„,0;3,0|re-'/'^'"^'|0„,l;3,0). (42) 

The precise form of the amplitudes depends on the choice of the proflle functions for the time duration 
of the interaction in the detector. In Ref. [72] two cases have been considered, (i) a proflle function 
for the detection narrowly centered around a detection time tn and (ii) a profile function with support 
on the time interval [^1,^2]) where ti is later than the time of emission, and ^2 is eventually taken to 
infinity. In case (i) the final expression for the probability is computed by integrating the amplitude 
Eq. (I42p squared, with respect to Id and by introducing a suitable normalization factor. In case (ii) 
the final expression for the probability is obtained by evaluating the amplitude Eq. (jl2]) squared at the 
limit ^2 — ^ 00. The results in both cases recover the standard oscillation formula. 

Regarding case (i) above, we note that the detection time to is also introduced as an external param- 
eter to the system and not as a variable, so our previous critique about the problems in integrating the 
squared amplitudes over time stands as in Sec. 4.2.3. Case (ii) is more unphysical, as the authors of Ref. 
|72j state, because it involves the limit t — )• 00. This case also corresponds to the treatment of neutrino 
production and detection in Ref. [75], even though in the latter reference the initial and final states are 
not discrete. The reason is that these works employ the first-order term in time-dependent perturbation 
theory, while at the long-time limit there exist significant higher-orders contributions to the total prob- 
ability. The leading-order terms dominate only at early times, when time-of-arrival interferences have 
not been built up. To see this, we note that this description of production and detection — i.e., treating 
the neutrino as an intermediate state — is formally similar to atom-field interactions in quantum optics, 
where photons cause transitions between energy levels in pairs of atoms. 
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In the quantum optics context, the vahdity of such approximations has been studied in detail |98j . 
and it is known that higher-order terms are necessary for the description of the long-time limit. The 
leading-order terms in the perturbative expansion of the finite-time propagator dominate only at times 
much smaller than emission time and their behavior cannot be extrapolated to the t — t- oo limit. For 
example, the Wigner-Weisskopf description of atomic decay [99], which is valid for long times, involves 
the solution of an integro-differential equation [98 1 , and it is equivalent to a resummation of an infinite 
number of Feynman diagrams from time-dependent perturbation theory |100j. 

5 The quantum measurement approach to particle oscillations 

In Sec. 4, we presented our review of the main problems in the existing theoretical approaches to particle 
oscillations. These problems originate from the lack of a first-principles derivation of the probability 
assignment relevant to the particle oscillation experiments. The probabilities defined in terms of squared 
amplitudes carry a time-dependence, which must be removed in order to obtain experimental predictions. 
The procedures employed so far for the removal of this time-dependence, such as the substitution of a 
'mean' time of arrival in the squared amplitude, or the integration of the squared amplitude over time, 
are mathematically unjustified. Moreover, they do not treat the time of arrival as a quantum variable. 
In this section, we present the quantum measurement approach to particle oscillations that aims 
to resolve the problems of the existing methods. The key point in this method is the construction of 
the probabilities relevant to particle oscillation experiments from first principles in quantum theory. To 
this end, we develop a general framework for measurements (Sec. 5.1) in which the time of detection 
(coinciding in the particle-oscillations context with the time of arrival) is a quantum variable, rather 
than a pre-determined parameter as in the common case of von Neumann measurements. This is a 
proper framework for the description of particle oscillation experiments, because the experimentally 
fixed quantity in the latter is not the time of detection but the location of the detector. The result is a 
non-standard expression for the particle-oscillation wavelength (Sec. 5.2). 

5.1 A general framework for quantum measurements 

In Sec. 5.1, we develop a general formalism for quantum measurements (to be applied in particle 
oscillations in Sec. 5.2) that treats the time of detection as a quantum variable. To this end, we generalize 
a method developed in Ref. [9], for the construction of probabilities starting with amplitudes that are 
defined at different moments of time. This method was applied to various problems, such as probabilities 
for the time of arrival, tunneling time |101] and for the study of non-exponential decays |102j . It contains 
ideas from the decoherent histories approach to quantum mechanics [29l |30l EU [32l |33l I103t I104J and it 
has many similarities to the Srinivas-Davies photo-detection theory |105) . 

5.1.1 Detection time as a variable in probability amplitudes 

The first step in our study is the derivation of a general formula for the probability of detection outcomes 
that does not treat detection time as a predetermined parameter. 

Let Ti be the Hilbert space of the quantum system under study; Ti may include the degrees of freedom 
of the apparatus in addition to the ones of the microscopic system. In the most general case, a quantum 
event, such as a particle detection, can be described as a transition of the system from one subspace of 
Ti to another subspace orthogonal to it. For example, 

• the emission of a photon from an atom corresponds to a transition from the one-dimensional 
subspace, defined by the electromagnetic field vacuum, to the subspace of single-photon states 
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• a von Neumann measurement corresponds to a transition from the subspace, in which the pointer 
variable X takes its pre-measurement values, to a subspace corresponding to possible measurement 
outcomes; 

• a time-of- arrival measurement through a detector, located at x = 0, may be idealized as a transition 
from the subspace corresponding to x G (— oo,0] to the subspace corresponding to x G [0,oo) f9]; 

• any particle reaction can be described as a transition from the subspace of states for the initial 
particles to the subspace of the product particles. In particular, a reaction in which a neutrino is 
annihilated can be described as a transition from the subspace of neutrino single-particle states to 
(other) lepton single-particle states. 

Hence, we assume that T-i splits into two subspaces T-Lj^ © 11- , where T-Lj^ describes the accessible 
states of the system, given that a specific event is realized, and Ti.^ is the complement of T-Lj^. We denote 
by P the projection operator onto %+ and by Q the projector onto %-. 

Once the transition has taken place, it is possible to measure the values of various observables 
through their correlation to a pointer variable. Let us denote by P\ projection operators corresponding 
to different values A of an observable that can be measured only if a detection event has occurred. For 
example, in neutrino oscillations. Pa niay correspond to a coarse-grained position variable for one of the 
product particles. The set of projectors P\ is exclusive {PxPx' = 0, if A / A') and exhaustive, provided 
a detection has occurred; i.e., J2\ P\ = P- 

In general, it suffices that the operators Pa be positive rather than projectors; they still have to be 
exhaustive on %+, but they do not have to be exclusive. 

We assume that the system is initially prepared at a state iV'o) G 'H+-, and that the dynamics is gov- 
erned by the self-adjoint Hamiltonian operator H. Our aim is to construct the amplitude iV'o! -^i [ii-, ^2]) 
that, given the initial state iV'o)) at some instant in the time interval [ti,t2]i a transition occurs with a 
recorded value A for the measured observable. We assume that [ti, t2] C [0, T], where T is a final moment 
of time, corresponding to a long integration time of an experiment. 

We first consider the case that the relevant time interval is small, i.e., we set ti = t and t2 = t + 6t, 
and we keep only terms of leading to 6t. The construction proceeds as follows. 

Since the transition took place within the interval [t, t + 6t], at times less than t the state lay within 
Ti-. This is taken into account by evolving \iPq) with the restricted propagator [106] into Ti-, i.e., with 
the one-parameter family of operators 

St = lim (Qe-^^*/^Q)^. (43) 

By assumption, the transition took place some time within the time interval [t, t + 6t] and at the 
end of this interval the outcome A has been recorded. This means that in the time-interval [t, t + 5t] the 
amplitude transforms under full unitary operator for time evolution e~*^ ~ 1 — i6tH. At time t + 6t 
the event corresponding to Pa is recorded, so the amplitude is transformed by the action of Pa, or of 

Y Pa) if -Pa is not a projector. Then, since there is no constraint, so the amplitude evolves unitarily until 
time T. 

At the limit of small 6t, the successive operations above yield 

|V'o;A,[t,t + <5t]) = -i6t e-'^^'^-''^PxHSt\^o). (44) 

We must emphasize here the important physical distinction on the role of time that is manifested 
in the derivation of the amplitude Eq. ()44p . The time of detection t is not identical to the evolution 
parameter of Schrodinger' s equation. Instead, it is a dynamical variable that determines the moment 
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that a physical event has taken place [7j. The construction of the amplitude \ip] X,[t,t + 6t]) above 
highlights this distinction: the detection time t is distinct from the time T at which the amplitude is 
evaluated. 

The amplitude {ipo'jX, [t, t + 5t]) is proportional to 6t, hence it defines a density with respect to time 
IV'o; A,t) := lim^t^o JilV'o; A, [t,t + 6t]), which by Eq. (04]) equals 

|Vo;A,t) = -ie-'^^^-'^PxHSt\^Po) = -ie~'^^Cx,tm, (45) 

in terms of the class operator 

Cx^t = e'^'^PxHSt. (46) 

Since the amplitude |V'o; A, t) is a density with respect to the time of detection t, integration over time 
is well-defined. In particular, this means that alternatives labeled by different values of t are disjoint 
|108j . Hence, we can integrate over t |107j . in order to obtain the total amplitude for detection happening 
at any time in the interval [ti,t2] as 



|V^;A,[ti,t2]) = -ie-^^^ f' dtCxMo)- (47) 



We note that if [P,H] = 0, i.e., if the Hamiltonian evolution preserves the subspaces T-L±, then 
IV'oJ A,i) = 0. When the Hamiltonian is of the form H = Hq + Hj, where [Hq,P] = 0, and Hj a small 
interaction, then to leading order in the perturbation 

Cx,t = e'^''^PxHie-'^°\ (48) 

Finally, we note that the amplitude Eq. (j44p has a more intuitive interpretation in the sum-over- 
histories formulation of quantum theory. If the operators P, Q and Pa are spectral projectors of con- 
figuration space variables, the amplitude Eq. (j44p corresponds to a path integral with a domain of 
integration restricted over a suitably defined class of paths [331 HlOt Hill I112J . Such path-integration 
domains are said to define spacetime coarse- graining s. 

5.1.2 Probabilities for detection time 

From Eq. (|47|) . we find that the probability p{X) that at some time in the interval [^1,^2] a detection 
with outcome A occurred equals 

p(A, [ti,i2]) = (V'; A, [tut2mX, [ti,t2]) = r dt f dt' Trie'^^'~''^PxHSlpoSt'HPx), (49) 

where po = \tpo){tpo\- 

Since the probabilities Eq. (j49p are obtained by the square of an amplitude, they are in general 
characterized by interference, i.e., the probability for an interval [tijts] = [^1,^2] U [^2)^3] equals 



p{X, [ii,t3]) =P(A, [tut2]) +p(A, [t2,t3]) + 2Re 



/ dt / dt'Tr(Cx,tPoCit, 

Jti Jt2 ^ 



(50) 



Probabilities are a measure of event counts; hence, the probabilities Eq. (j49p cannot describe the 
outcomes an experiment, unless the terms 



2Re 



dt f ' dt'Tr (Cx,tPoClt, 



t2 rts 
ti 
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vanish. A measurement apparatus is a macroscopic system with a large number of degrees of freedom, 
therefore one expects that consistent probabiHties can be defined given a substantial degree of coarse- 
graining J3HI32]. In particular, there is a coarse-graining time-scale r, such that the non-additive terms 
in Eq. ([501) ^^^ strongly suppressed if |f2 — ti\ >> r and |i3 — t2\ » t. Then, Eq. (|19|) does define a 
probability measure when restricted to intervals of size larger than r. We explain how this coarse-graining 
scale arises in the Appendix A, where we consider a model measurement scheme. 

Note that we are interested in constructing a POVM that describes measurement outcomes, i.e., the 
physical system under consideration includes the measuring apparatus, which is a macroscopic body 
with a large number of degrees of freedom. The pointer variables are essentially classical and classical 
fluctuations (thermal or statistical) impose a lower limit r to the temporal resolution for the system. 
Hence, there is no meaning in discussing about the probabilities associated to time-intervals of width 
smaller than r (see the discussion in Sec. 3.2). 

Temporal smearing. Assuming a finite resolution scale r, one may use the amplitudes Eq. ()47p to 
define a POVM. It is more convenient to smear the amplitudes at a time-scale defined by r rather than 
considering sharply defined intervals as in Eq. ([49]). 

To this end, we introduce a family of functions /t(s), localized around ,5 = with width r, normalized 
so that limT-_s.o fris) = 26{s) |109j . For example, we may employ the Gaussians 

We then define the amplitude \ilJo;X,t)r, localized around the value t with width r, as 

IV'O; A, t)r = I ds^fris - t)|V'0; A, s) = Q^A.r I V'o) , (52) 

where we wrote 



Ct,x,r = j ds^fris- t)Cx,s. (53) 

The square amplitudes 

Pr{X,t)=r{^o;X,t\i;o;\t)r = Tr[cl^JoCx,t,r) (54) 

are of the form Tr[/)onr(A, t)], where 

UriX,t) = Cx,t,rCl,,. (55) 



The positive operator 



n^(Ar) = 1 - / dt I dXUr{X,t), (56) 



corresponds to the alternative M that no detection took place in the time interval [0, oo). nT-(A^) together 
with the positive operators Eq. ([55p define a POVM on ([0,oo) x 17) U {A/"}, where O is the space of 
possible values of A. The POVM Eq. (I55p determines the probability densities that a transition took 
place at time t, and that the outcome A for the value of an observable has been recorded. 

In Eq. (155p the time of arrival is a genuine observable (a random variable) and not a parameter, so 
integration over time in Eq. ([55]) is well defined and it leads to a partial POVM IIt-(A) = J^ dt IIt-(A, t). 
This equals 

/•oo roo 

flr{X)= ds dsGi{s',s)Cx,sCi,„ (57) 

JO Jo 
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in terms of the kernel 

Gl{s', S) = P dt^fris - t)fr{s' - t). (58) 

If the POVM acts on states with support on values for the time-of-arrival much larger than r, the 
Gaussian smearing functions Eq. (j5ip can be used. In this case, 

(«-s')2 
Gi{s,s')=e i^^. (59) 

Finally, we note that if the temporal resolution r of the detection is much larger than any time-scales 
in the amplitudes Cx^tlipo), then smearing is not necessary. Then, simply consider the probabilities Eq. 
I9|) for the interval [0,oo), to obtain 

p(A) = r dt r dt' Tr (Ci ,poCx,t') . (60) 



Appendix A contains an application of the above formalism to a system consisting of a microscopic 
particle and a measurement apparatus. Then, we show that the apparatus defines a time-scale Tdec^ such 
that a choice of r in the smearing functions of the order of Tdec, guarantees that the probabilities in Eqs. 
()55p and (j57p are well defined. Hence, Tdec sets a scale of minimal temporal coarse-graining necessary 
for the quasi-classical attribution of definite macroscopic properties to the pointer variables. Thus it 
corresponds to the maximal temporal resolution of the detector. In the Appendix A, we also show that 
von Neumann measurements are a special case of the framework we developed in Sec. 5.1. 

5.2 Particle oscillations 

5.2.1 The construction of a POVM for particle oscillations 

We now apply the formalism developed in Sec. 5.1 to the treatment of particle oscillations. The 
oscillating particles are detected by means of their reactions, hence, it is necessary to include into the 
description the degrees of freedom of all particles participating in such processes and to employ quantum 
field theory. 

In what follows, we will not specify the type of oscillating particles, so our results are valid for both 
neutrino and neutral boson oscillations. To this end, let us denote the oscillating particles as A and 
assume that they are detected by means of the process 

A + Bi + ...Bm^Di... + Dn, (61) 

where Bm, Dn are particles (different from the A particles), labeled by the indices m = 1, . . . ,M and 
n = 1, . . . , A^. In this description, we assume that the A particles are produced in a production region, 
they propagate and they interact with (one or more of) the Bm particles, which are located in a detector 
away from the source. The interaction of the A particles with the Bm particles produces the products 
particles Dn, which are the ones that are being detected. Relevant observables are the Dn-particles' time 
of detection, position, momentum and so on. These observables are determined through macroscopic 
pointer variables in the detector. 



The Hilbert space. The Hilbert space Titot, in which the process Eq. (I6ip are described, is expressed 
as a tensor product Titot = T~{-a (X" Tir- Ti-A is the Fock space F{Hia) corresponding to the A-particles. 
Explicitly, we write 

-F(Hia) = C®HiA® (HiA Hia)s,a e {HiA ® Hia ^ Hia)s,a . . . , (62) 

30 



where S refers to symmetrization (bosons) and A to antisymmetrization (fermions) and Hu^ is the 
Hilbert space describing a single A particle. The single-particle Hilbert space is a direct sum (BiTii of 
the mass-eigenspaces Tii . The degrees of freedom of the Bfn and Dn particles are incorporated into the 
Hilbert space Tir- In general, Tir is a tensor product of Fock spaces, one for each field other than A, 
participating in the process Eq. (|6TI) . 



Events. In order to apply the formalism of Sec. 5.1, it is necessary to identify the subspaces T-L± that 
define the transition under consideration. Since now the detection proceeds through the process Eq. 
(j6ip and the A particles are not directly observable, the transition corresponds to subspaces of Tir- The 
Hilbert space Tir is decomposed as Tio © Tiprod, where ^o is the subspace of states prior to the decay 
A + Bjn — >■ Dn, and Hprod is the subspace corresponding to states of the decay products. This means 
that we identify Ti- with "Ho and 71+ with Tiprod- 

The following are examples of such subspaces for specific detection schemes. 

• If the Kl boson is detected by means of the decay K^ -^ Svr*^ , 'Hr is identified with the Fock space 
Jvg for vr" particles, Tio is the vacuum subspace of J^g and Tiprod is the subspace of J>g with 
non-zero particle number. 

• For the detection of neutrinos through the reaction Uf, + p — t- n + e~ , we have only zero- and 
one-particle states for each of the p, n, e~ particles. Therefore, the relevant Hilbert space for each 
particle p, n and e~ is the direct sum of the corresponding vacuum subspace T-L^p^^ and the single- 
particle Hilbert space 'Hn^p^e- The Hilbert space Tir is then the tensor product {'Hp'^'^ ® Tip) (8> 
{'Hn'^'^ © T-Ln) © {Ti'^'^'^ © Til)- If the detector records product electrons, the subspace Tiprod equals 
{Til"" © Til,) ^ ^'^'"n^ ® ^n) ® ^e and Tio = {Ti;"" © Ti^ © {Ti^^" © 7^^) © Ti^^. If both product 
particles are detected then Tiprod = iTip^'^ © Tip) © Ti\ © Ti].. 

• If the neutrino is detected through scattering, for example in a reaction v + e~ — )• v + e^ , then 
the Hilbert space Tir is identified with the single-particle Hilbert space Ti], of the electrons, which 
splits into two subspaces Tio®Tiprod- Tio corresponds to the spectral projector of the Hamiltonian 
in the range [0,So]) and Tiprod to the spectral projector of the Hamiltonian in the range [£'o,oo). 
£^0 is the maximum energy in the initial state of the electrons. 

We will denote as P the projector P : Tir — ^ Tiprod and Q = 1 — P its complement. The corresponding 
projectors on Titot are 1 © P and 1 © (j. 

Measurements carried out on the product particles correspond to a family of positive operators 1©IIa 
onto Tiprod, where Ha are positive operators on Tir corresponding to different values A of the measured 
observables. They satisfy the completeness relation / dXIl\ = P. 

Dynamics. We assume a Hamiltonian of the form 

H = HA(^i + l(^Hr + Hi, (63) 

where Ha is the Hamiltonian for free A particles, Hr is the Hamiltonian for the Bm and Dn particles, 
and Hj is the interaction term. 

In the following, we will consider a single-particle initial state for the A particles, so that the Hamil- 
tonian Ha is restricted on the TiiA subspace, i.e.. 



i^A = Ev"^? + p'^^' (64) 
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where p is the ^-particle momentum, and Pi the projector onto the eigenspace corresponding to particle 
mass rrii. For ease of notation, we ignore the A-particles' spin degrees of freedom since these do not 
significantly contribute the results pertaining to oscillation probability. 

We assume that any particles B^ that are present prior to detection are almost stationary. This 
condition defines our reference frame, that in most cases coincides with the laboratory frame. Assuming 
that the initial state for the Bm particles has support to values of momentum much smaller than their 
rest mass, the restriction of Hr in the subspace TIq is a constant. It is convenient to set this constant 
equal to zero. Then, 



Hr 



€th + J2^ 



where 



Mh„ + Pl 



MdJ 



etfc = ^MB„-^MB„ 



P, 



(65) 



(66) 



is the threshold energy of the A particles for the process A + B^ -^ Dn- In the above, Mb^ and Md^ 
are the masses of the particles Bm and Dn respectively and p^ are the momentum operators for the Dn 
particles. 

We specialize to the case that the ^-particle is annihilated in the reaction through which it is detected. 
We will therefore consider an interaction Hamiltonian of the form 



Hi 



E 



d^x 



■6,(x)^„,J+(x) + 6|(x)[/^J-(x 



(67) 



where bi,bl are annihilation and creator operators on T-La, i labels mass eigenstates, J^ (x) are current 
operators of flavor a defined on Tir, and Uai is the mixing matrix. A more general case (relevant 
to bosons) would involve a kernel instead of a constant for the mixing matrix Uai, reflecting the fact 
that the mixing coefficients may depend on momentum. Here, however, we shall consider initial states 
sharply concentrated in momentum, for which a constant value of Uai provides a good leading-order 
approximation. The case of detection through scattering involves a different interaction Hamiltonian 
quadratic to the A-field variables and it will be treated in a different publication. 

The current operator J^ involves products of annihilation operators for the B particles and creation 
operators for the D particles. Since no ^-particles are created during the detection process, the initial 
state |0o) in Tir must satisfy 



4(x)|0o) 



0. 



Amplitudes. We assume an initial factorized state for the total system \iPq) 
T-La is a single-particle state 



IV'o) = E/^'^^1WV'^o(x)|o)a, 



The state |V'o) on 



(69) 



where \Q)a is the vacuum of the Fock space T-La- We only examine single-particle states for the A particle 
in this paper, so the choice of fermionic or bosonic statistics makes no difference. The state |0o) is, in 
general, a many-particle state, since the detector contains a large number of i?-particles and one cannot 
specify a priori which one participates in the reaction. 
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We now construct the amplitude {ipQ, (J)q; a, X,[t,t + St]) , for the detection of a particle produced from 
a reaction through flavor a, at the time-interval [t, t + 6t] and with a measurement outcome A. From 
Eq. (flSl) . we obtain 



\iPo,cl)o;a,Xdt:t + St]) = -i6te-'^^Ca,xMo) (^\M, (70) 

in terms of a class operator Ca,\,t- 

Since [Ha ®1 + 1® HrA®Q] =0i we can use Eq. (Ii8|l in order to obtain the amplitude to leading 
order in perturbation theory. We find 



i 

where 'ilji{x,t) is the solution of Schrodinger's equation for the ^-particle Hamiltonian, with initial 
condition tpiixjt) = ipio{x). 

Probabilities. We next employ the amplitudes Eq. (j7ip for the construction of a probability density, 
using the methodology described in Sec. 5.1. To this purpose, we specify the temporal resolution Tdec 
of the detection scheme and we employ Eq. (j55p . Oscillations are observed in relation to the distance of 
the detector from the source and not directly with respect to the values of time of arrival — whether the 
time of arrival can be distinguished or not. Hence, we consider the partial POVM Eq. (j57p . obtained 
by integrating over the time of arrival in Eq. (|55p . 

The probability that a decay through flavor a has happened at some time in [0, T] and that the value 
A for an observable of the product particles has been found equals 

p^{X) = J2 f dt I dt' I d^xd^x'i)*{-s:,t')ijj{-iL,t)U*^iUajGi{t' - t)Rl{x,^,t- t'), (72) 



ii^(x,x',t-tO = (c/>o|J-(x')VnIe-*^'-(*'-*)VnIj+(x)|</.o). (73) 



where 



and Gi{s' — s) is the kernel Eq. ([55]) . for r = Tdec, that carries the dependence on the resolution scale 
Tdec of the detector. 

Observables. In order to complete the description of the measurement scheme, we must specify the 
observables that are being recorded at the detection of the decay particles. This is encoded into the 
specification of the positive operators IIa in Eq. ([73]) . In the most general case, the positive operators 
Ilx correspond to joint measurements of phase-space variables, i.e., the apparatus records the position 
and momentum with an accuracy corresponding to a phase space volume of size much larger than h ' 
|29l I115J . However, the degree of coarse-graining associated to the sampling of momentum is much 
larger than the coherences of momentum in the quantum states, relevant to particle oscillations, so for 
simplicity we specialize here to the case of sampling only with respect to position. In Sec. 5.4.2, we 
show that this choice does not affect the main conclusions. 

We assume that the positive operators H^ correspond to the sampling of the position X of one of the 
product particles, with an accuracy of order 6. For example, we can use Gaussian operators of the form 

nx= /d3X'e"^^i^|X')(X'|®l, (74) 
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where the tensor product with unity refers to the remaining degrees of freedom in %prod- 

The vector J^{^)\(J)q) in Eq. ([73]) corresponds to the state of the product particles for decays that 
have taken place in a neighborhood of the point x. In order for the operators IIx to provide a definite 
determination of the locus of the decay event with an accuracy 5, it is necessary that they have J^(x)|(^o) 
as approximate eigenstates, i.e., 

fl(X)J+(x)|(/)o) ~0,if |x-X| »5 
n(X) J+(x)|0o) - J+(x)|0o),if |x - X| « 5 (75) 

Heuristically, this means that the width 5 of the sampling is much larger than any length scale char- 
acterizing the interaction. As we discussed in Sec. 3.1, the accuracy (5 of a position measurement is 
macroscopic, because a definite recording of position is only possible at a scale coarser than the fluctu- 
ations of the corresponding pointer variable. On the other hand, the currents -/^^(x) reflect microscopic 
(in fact sub-atomic) length-scales. Hence, there is a sharp separation of scales and the conditions Eq. 
()75p are expected to hold. Moreover, at this level of coarse-graining, the localization of a single prod- 
uct particle determines the location of all other product particles at the moment of detection; their 
wave-functions after the decay have the same support in position as J^(x)|(/)o). 

While 5 is much larger than any microscopic scale characterizing the decay event, a particle oscil- 
lation experiment is only meaningful if 6 is much smaller than the length scales of particle oscillations. 
Therefore, to a first approximation, we can substitute ■i/'a(x, t) by ■0a(X, t) in Eq. (172]) . 

There are three length scales entering the expression for the particle oscillation probability: a mi- 
croscopic scale Imic associated to the decay process, a macroscopic scale 5 associated to the detection, 
and a macroscopic scale L characterizing the distance of the detector from the source. These scales 
satisfy the condition l^ic « ^ « L, that results into the kernel R being approximately proportional 
to (5(x, X)5(x',X). In order to estimate the factor of proportionality and its time-dependence we note 
that, by virtue of Eq. (175p . the kernel R becomes 

iZ^(x,x',t-t')^(0ol4-(x')e-^^^(*'-*)j+(x)|<^o), (76) 

where x and x' are assumed to lie within a distance of order 5 from X. The vectors J^{x.)\(po are then 
proportional to |x)i (g) |x)2 (8> . . . fS" |x)„, where by |x)j we mean a state for the product particle labeled 
by n, localized around the point x. At a scale that we can approximate -Rx by a local kernel, the matrix 
elements are given by 



5o 



J^(x')e-^^'-(*'-*) J+(x)|0o) ^ i^ie^^*'^(*-*') X n f / dx f dx' (x|e'^v/^^+^(*'-*V)) , (77) 



where Ki is a constant, Cx,& is a region around X of width 5, and (x|e~ v ^^ c„v - '\-x') is the 
propagator for a relativistic particle of mass Mq^. Substituting the integration over Cx,<5 with an 
integral with a Gaussian smearing function, we obtain the local approximation to the kernel 

i?i(x,x',t-t')- i^i<5^(X,x)5=^(X,x')e"*'^("')G2(t'-t). (78) 

The function G2{s) arises from the smearing of the free particle propagators, in a region of width 5 

G2{s) = \{ld^P e-(V^^C^-A^^n)-5^p^ (79) 
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The index n runs over all product particles D^- If ^£,^5 » 1, the saddle-point approximation 
applies and 

^-'-'-n ( ,,.(, ^,,,.^,J '^ (80) 

The function G2{s' — s) suppresses contributions from the amplitudes tpa{^,s) at different times s 
and s' only if |s — s'| > Tgup, where the suppression time scale Tsup = MminS'^ is determined by the 
mass Mmin of the lightest product particle. The lowest conceivable value of Tgup is obtained by setting 
Mmin = f^e and taking for 5 the atomic time-scale Q.lnm. We then obtain Tsup > W~^^s. A more 
realistic value for r^up is obtained by placing 5 in the macroscopic regime, since it corresponds to the 
macroscopic accuracy in the determination of a particle position. For 6 ~ 100n?7i, we find Tgup ~ 10ns 

[m]. 

In the regime relevant for particle oscillations, the precise form of the kernel 6*2(5) does not affect 
significantly the particle detection probability. The most important information contained in ^2(5) is 
the suppression time-scale Tg^p than is uniquely determined on dimensional grounds. Hence, none of the 
approximations used in the derivation of Eq. (I80p affects the main results. 

Eqs. ()72p and ()78p define a probability density for the particle's position x, at scales much larger 
than 6, 

Pa{x) = K [ dt [ dt'Aa{t,x)Alit',x)e-'^^*'-^'^F{t',t) (81) 

JO Jo 

where Aa{t,x) = J2i Uaitpi(t,x) is the probability amplitude corresponding to the flavor a, and 

F{t' -t) = Gi{t' -t)G2{t' -t). (82) 

The probability in Eq. ([8T]) can be explicitly computed for any type of oscillating particle provided we 
specify an initial state. 

5.2.2 Wavelength of particle oscillations 

Next we specify a class of initial states for the A particles relevant to the description of oscillation 
experiments. It is convenient to assume that the solid angle connecting the production region to the 
detector is very small so that only particles with momentum along the axis that connects the two regions 
are detected, and the problem becomes one-dimensional. We employ a /3-flavor superposition of Gaussian 
mass-eigenstates, centered around x = 0, 

i^o{-) = T.Uhj^.-'^^'''^^ (83) 

where a is the spread of the Gaussian that we take to be a free parameter. The dependence on the 
mixing matrix Ua- is a consequence of the creation of the -B-particles through an a-flavor current. The 
Gaussian form of the initial state is chosen for convenience. In Sec. 5.4.2, we will show that the results 
obtained from a state of the form Eq. ()83p hold for a much larger class of initial states. 
For the initial state Eq. ([831) the amplitudes Aait,x) in Eq. (fSTI) are 



dp _f2 
— f 
27r 



Aait,x)=Y,U*piUaii'iTTa^)^/^ / ^g-^(p-p,)2+ipx-^i?,(p)i-r,(p)t^ (84) 
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where Ei{p) = \/mf + p^ and Ti{p) are the decay rates in the different mass eigenspaces |117| . 

We calculate the amplitude in the regime where wave-packet dispersion is negligible. To this end, 
we employ the approximation of expanding Ei (p) to first order in p — pi and F j (p) to zero-th order in 
p-Pi, i.e., 

E^{p) = Ei + v^{p - Pi); Ti{p) = Ti, (85) 

where Ei = Ei{pi), Vi = pi/Ei, and Tj = Ti{pi). 
We then obtain 

Substituting Eq. (f86|l in Eq. ([8T|) and letting T — )• oo, we obtain the probability density for the detection 
of flavor a at X = L 

p^{V) = K I ^5,e"'^^ + ^2e"^^"'^^''i?e(I-,- e^'^-^) 1 , (87) 



where 



i<j 



12 ^ ^~^ r .l>:P(^\.-^i'~i(E^-^tH)i 
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and 



xe — 5^2-^ -2i^»-^iK-H^ fthK^ (89) 



h, = ^^^^-^^^^-{m-P,). (90) 



In the derivation of Eqs. (j87H90p we have used the fact that L » a, that allows us to take the 
range of time integration into the whole real axis. 

The expressions Eq. (j88p and Eq. (|89p simplify significantly in the physically relevant regime where 
the time-scale of the oscillations is much smaller than the characteristic timescales Tsup and Tdec- We saw 
that the physically relevant values of Tsup is of the order of nanoseconds, while the resolution time of a 
solid state device is also of the order of nanoseconds. Hence, the timescales corresponding to Eq. (j90p 
are many orders of magnitude smaller, either for neutrinos or for neutral bosons. This condition also 
implies that the energy and decay constant differences \Ei — Ej\, \vi — Vj\ and \Ti — Tj\ are very small 
in comparison to Ei, Vi and Fj, respectively . We can therefore take their values, in all non-oscillating 
terms of Eq. ([88]) and Eq. ([89]) . as approximately equal: Ei c^ E, Vi = v and Fj ~ F. We then obtain 



Si = |C/*^f/„,|2^_^F(0)e ^ (91) 



V 



T., = C/*,C/,,C/^^.[/*^.il!^^^F(0)e ^ e^^ ^ (92) 
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The oscillating term at the end of Eq. ()89p does not depend on L and the corresponding phase is 
much smaller than kijL (it vanishes for neutrinos) so we can ignore it. Then, Eq. (j87p becomes 



p„(L) = K'y£ \U*p,U^,\h-^^'' + J2 2^e lu*piUaiUpjU*je ^"^^^^ ^^e*^^) j , (93) 

where terms in Eqs. (|91ll92p have been absorbed into a redefined constant K'. Eq. (j93p is similar in form 
to the common expression for particle oscillations Eq. (j20p and it allows a direct comparison: it predicts 
a different value for the oscillation wavelength. For example, if eth/Ei << 1 then kij = mf/pi — m'^/pj, 
i.e., the value of kij is twice that of the standard expression. 

« p, where p is the mean momentum, Eq. 



In general, assuming that \pi - 


- pj\ << p and \mi — mj\ 


pi]) yields 






h, = {l-'^){m^-m] 



(94) 



In the ultra-relativistic limit Eq. (|94p yields 



kij = il-^){m^-m])/E. (95) 

The oscillation wavelength carries a strong dependence on the threshold energy; this dependence does 
not arise at all in the standard treatment. 

It is important to emphasize the robustness of this result. While the first-principles treatment leading 
to Eq. (|93p is involved Eq. (j93p itself can be recovered by simple heuristic arguments. The key points are 
(i) that the summation over different times takes place at the level of amplitudes, (ii) that the time-scales 
characterizing the oscillation are much smaller than the characteristic time-scales of the measurement 
scheme, and (iii) that the presence of an energy threshold involved in the detection is equivalent to the 
presence of a constant potential V = —eth, in which the j4-particle propagates prior to detection. With 
these three conditions, Eq. (|90p can be derived even with the plane wave method (Sec. 4.1.1). 

Separation of timescales. The key observation leading to this result is the large separation of the 
time-scales relevant to the problem. The terms responsible for the difference between the present result 
and the usual derivations arise from the contribution of the interference between different times of arrival 
in Eq. ()8ip . These take place at a time-scale of order 6tij ~ | — — -^ | . In the ultra-relativistic regime 

Am?. 

relevant to neutrinos, 5tij ~ L „-p , i.e., 6tij is a small fraction of the classical arrival time td — L. For 
L ~ Ikm, E ~ lOOMey, and Am?- ~ leV , we find 5tij ~ lO^^^s. On the other hand, the characteristic 
suppression timescale is Tsup ~ rne^"^, that is of the order of 10~®s, and the temporal resolution of the 
detector T^ec is of the order of 10~^s. There is a sharp separation in the orders of magnitude between the 
decoherence time-scale and the interference time-scale, which means that time-of-arrival interferences 
are not suppressed and their contribution to the oscillation wavelength persists. 

5.3 Critical analysis of the quantum measurement approach to particle oscillations 

We proceed to an analysis of the method presented in Sees. 5.1 and 5.2 in order to show that it resolves 
the problems in the theoretical treatments of particle oscillations as presented in Sec. 4 and that it does 
not suffer from the problems associated to derivations of non-standard formulas for particle oscillations. 
We will also discuss the points at which the method presented here needs further improvement. 
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The quantum measurement approach to particle oscillations was constructed with the avowed aims 
of providing a derivation of probabilities relevant to particle oscillation experiments that is based on the 
first principles of quantum theory. It also incorporates a fully quantum treatment of the time of arrival 
as a physical observable. To this end, we go beyond the modeling of experiments by von Neumann 
measurements that are inappropriate for the treatment of particle oscillations. 

The above aims are fulfilled by the general framework for quantum measurements developed in Sec. 
5.1. In this framework, detection time is treated as a physical observable thus incorporating the crucial 
distinction between time as a parameter of Schrodinger's equation and time of arrival. The end result is 
a class of POVMs that are applicable to a large class of measurement situations of which von Neumann 
measurements are a special case. 

In Sec. 5.1, we present our method in its most general and abstract context. The method applies to 
any physical system and the only inputs are the Hamiltonian operator and the specification of observables 
that are being measured. In Sec. 5.2, the method is applied to the special case of particle oscillations. 
The result is a probability distribution pa{t,L) for detection events depending on time of arrival t and 
distance L. Since the time of arrival t is treated as an observable, integration of the probabilities over 
t is well-defined and unambiguous. It leads to a probability density Eq. (|93l) for detection, solely as a 
function of the distance L. 

The probability density Eq. ()93p is constructed through a POVM, hence, it is a linear functional of 
the initial state po. The construction thus respects the convexity properties of quantum states and the 
resulting probabilities have no ambiguity in their statistical interpretation. In particular, the constant 
K' that appears in Eq. (I93p is not inserted by hand for normalization purposes, but it appears as a result 
of an evaluation of the kernel Eq. ()73p . Its precise value can be fully computed once explicit expressions 
for the positive operators IIa and the initial state have been provided. 

In the physically relevant regime of parameters, the resulting probabilities take a particularly simple 
form that could also be obtained by other means. For example, an equation analogous to Eq. (I93p 
can be obtained by the plane wave method using an unequal-times prescription. The key difference 
is that our result is a last step of a precise procedure, starting from first principles and making no 
assumptions unwarranted by quantum theory. Hence, even if our expression for oscillation wavelength 
(partially) coincides with other non-standard oscillation formulas, the methodology and the underlying 
assumptions are completely different. In particular, our derivation requires no assumption about the 
relation of energies Ei and momenta pi, in the different mass eigenspaces of the initial state. In our 
opinion, such assumptions about the initial state require an artificial fine-tuning of parameters that 
cannot be justified by the physics of the production process. A reliable physical prediction should be 
insensitive to tiny details in modeling of the initial state. Our results satisfy this property: we show in 
Sec. 5.4.2 that Eq. ()93p is valid not only for initial states of the form Eq. ()83p but for a large class of 
initial states sharing some broad qualitative characteristics. 

The use of the unequal-times prescription led to a prediction of oscillations for recoil particles, for 
instance the muon in vr — t- /iz^ or the A in 7r~p — t- AK^. Here, in contrast, we defined the POVM 
for particle oscillations in terms of the detection of product particles, since the neutrino is not directly 
observed. Hence, we can apply Eq. (f93|) also for the detection of the A particle. This has a single mass 
eigenspace, therefore no off-diagonal terms and, trivially, no oscillations appear. 

Eq. (j90p for the oscillation wavenumber followed from a lengthy calculation from first principles 
of the quantum formalism. Nonetheless, its difference from the standard oscillation formula can be 
simply explained: the standard result involves an assumption that the contribution of time-of-arrival 
interferences to the total probability is suppressed. This assumption is implemented in different ways 
in the literature, but it underlies all existing derivations. We do not make such an assumption here. A 
regime corresponding to suppressed time-of-arrival interferences also exists in our formalism but it does 
not correspond to the physical range of parameters. To see this consider Eq. ()8ip . If the kernel F{s — s') 
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Figure 1: The dependence of the oscillation wavelength A on energy E is shown in this plot, where the ratio 
XAm^/E is plotted as a function of E/eth- Ato^ is the neutrino squared masses difference and Cth the energy 
threshold for the detection process. The dashed line represents the standard expression, the dotted line the 
non-standard expression Eq. ([2]) and the solid line our result, Eq. (|95|) . 



can be substituted by a delta function the resulting probabilities are proportional to ^ dt\Aa{t,L)\^, 
and this expression leads to the standard result for the oscillation wavelength. However, this condition 
requires that either the decoherence time-scale Tdec or the suppression time-scale Tgup in F{s — s') is 
much smaller than 5tij, which is not the case for particle oscillation experiments. 

Different experimental prediction. The quantum measurement approach to particle oscillations 
makes a definite physical prediction, according to Eq. (|95p . the product kijE increases with energy. In 
particular, it varies by as much as 50% from its smallest value, for energies near the threshold ej/i, to its 
largest value, at the limit E » eth- In contrast, the standard expression predicts a constant value of 
kijE. Previous non-standard oscillation formulas also predict a constant value of kijE and they cannot 
be distinguished from the standard one, in absence of an independent determination of neutrino masses. 

The dependence of the oscillation wave- number kij on energy, in Eq. (j95p . involves a single unknown 
parameter, because the threshold energy for each detection process is a known quantity. Hence, Eq. (I95p 
can be directly compared with the standard expression in terms of best-fit to existing data. Moreover, 
if the accuracy of energy sampling cte is not larger than eth, then the corresponding distributions of 
number of neutrino detections as a function of energy, have a significant statistical distance. As such, 
the two distributions can be sharply distinguished given a sufficient number of events. 

One way to see this is by considering the two-neutrino case where the only unknown parameter is 
the difference Am^ of squared neutrino masses. We further simplify the description by ignoring the 
dependence on energy of the pre-factor K' on energy, so that the energy dependence of the number of 
events is dominated by the behavior of the oscillating term. 

Then, according to the standard oscillation formula, the number of events as a function of energy is 
given by a distribution 

n,{E)^ sin' ^, (96) 



while we predict a different distribution 



77-2 (-E) ~ sin 



2E ^ 2E' 



(97) 



Since Am^ is unknown, we define a parameter /x^, as the value of the product kE at the limit E — )• oo. 



Hence, fi' = ■^/!\m? in Eq. (j96p and ^^ = /\m? in Eq. (j97p . The event count at high energies determines 
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Figure 2: Heuristic plots of the distributions pi{x) Eq. ([M]) (dashed line) and P2{x) Eq. (|^ (sohd line) as a 

2 r 

function of x = E/cth, for different values of the parameter a = ^ — . The distributions are normalized to 1 for 
all events with x & [1, 10]. 




Figure 3: The statistical distance di2 = 8(1 — J dx^pi {x)p2 {x)) of the two distributions pi (a;) andp2{x) provides 
a measure of their distinguisha 
the statistical distance with a. 



T 

a measure of their distinguishability. Here it is plotted as a function of a = ^ — . There is an overall increase of 



a value of /i, and the predictions of Eqs. (f96|) and (f97|l can be compared at low energies (as E approaches 
eth) where they diverge. 

Defining a = ^^ and introducing the scaled energy x = E/eth, we compare the probability distri- 
butions 



Pi[x) 



Ci sin — , 

X 



P2{x) = C72sin2^iL_J-l. 



(98) 
(99) 



The constants Ci and C2 above have been introduced for the purpose of normahzation over recorded 
events in any specified energy range. 

The difference in the two distributions at low energies is highlighted in Figs. 2 and 3. Fig. 2 is a plot 
of the two distributions as a function of x for different values of a. Fig. 3 plots the statistical distance of 
the two distributions as a function of a. 

At low energies the divergence of the present oscillation formula from the standard one is stronger. 
This fact is rather significant, because it is at low energies that the discrepancy of the number of observed 
neutrino events has been observed, in the MiniBooNE experiment. The result obtained in this paper is 
a possible explanation of this discrepancy that does not require the introduction of physics beyond the 
Standard Model. 

Open issues in the formalism. The main advantages of the quantum measurement approach to 
particle oscillations are that it makes no assumptions extrinsic to quantum theory, and, by employing 
sound probabilistic reasoning, and leads to well-defined probabilities. 
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There are, however, two points where further improvement is needed. First, we need to guarantee 
the uniqueness of the results, namely, we must show that alternative methods for constructing joint 
probabilities for time of arrival and position would lead to the same conclusions when applied to particle 
oscillations. We believe that the simplicity and robustness of our final result is a strong indication that it 
would persist in alternative treatments. Nonetheless, several different approaches for the quantum treat- 
ment of the time of arrival exist, and their generalization to the particle-oscillation context could provide 
an independent theoretical confirmation. It is to be noted, however, that many of these approaches to the 
time of arrival — especially ones relying on probability currents — cannot guarantee positive probabilities 
when the time-of-arrival interferences are significant (see. Sec. 4.2.6). Hence, they may not be applicable 
to particle oscillations. In our opinion, the most promising candidates for an alternative treatment of 
the time of arrival are the axiomatic method of Kijowski ^23j and the decoherent histories approach of 
Halliwell and Yearsley |56] . 

An unsatisfactory point of our approach involves the probability of no-detection, Eq. ()56p . While 
it is mathematically well-defined, and its meaning is physically clear, i.e., not all neutrinos generated 
at the source will be detected, the definition Eq. (j56p follows from the normalization of the POVM 
and it is not obtained by independent arguments. The problem with an independent definition of this 
probability is that the natural candidate, namely, limt_^oo Tr^StpoSJ), is identically equal to 1 [9]. This 
fact is a special case of the quantum Zeno effect |113) that is generic in quantum theory. The Zeno effect 
creates problems for any direct definition of the probability of no-detection. The usual construction of 
the probability of no-detection proceeds through the computation of a state's persistence amplitude. 
This method is reliable for the study of exponential decays but it does not constitute a fundamental 
definition, because it cannot guarantee positive probabilities |114j . The definition Eq. (j56p for the 
probability of no detection has been shown to be physically sound in other applications of the formalism 
[Hi 1101] and it reproduces the standard result in cases where independent methods for its determination 
exist (exponential decays of unstable states |102j ). A better definition is desirable for the conceptual 
completeness of the method but the present one, Eq. (j56p . is mathematically and physically sound and 
it does not curtail the applicability of our formalism. 

5.4 Application of the formalism to more general settings 

Starting from Eq. ()55p . Eq. ()90p in Sec. 5.2 was derived using three approximations. 

• The local approximation Eq. (I78p for the kernel R in Eq. (I72p . 

• The choice Eq. (I83p for the initial state of the A-particles. 

• The evaluation of the amplitudes Aait,^), in the limit of negligible dispersion. 

We now proceed to show that our result, in particular, the expression for the oscillation wavelength 
is insensitive to the above approximations. 

5.4.1 More general samplings 

We first examine the case of samplings IIa on the detected particles, that are more general than the ones 
considered in Sec. 5.2.1. In the process, we will also elaborate on the validity of the local approximation 
Eq. (f78l) that was employed in the derivation of Eq. (i90|) . 

To this end, we return to Eq. ()72p and, instead of position samplings on one product particle, 
we assume that the particle's momentum is also measured. The positive operators P\ will therefore 
correspond to (coarse) joint measurements of phase space observables, i.e., the position X and momentum 
P of the recorded particle. 
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An important property of phase space POVMs is that they are practically equivalent at the limit 
of large coarse-graining |29^ I115| . This limit refers to properties of macroscopic observables, such as 
the pointer variables in a measurement scheme. In particular, if Jlu and 11^ are two different positive 
operators, corresponding to the same phase space region of volume [U] » Ti"'^, then 

where n is the dimensionality of the phase space |118j . 

We consider Eqs. (f72]l and ([73]) with A corresponding to the pair (X, P) and IIa := Hx.P) a POVM 
that corresponds to joint position and momentum measurements with respective accuracies 5 and Sp, 
such that Sp6 » h. We evaluate Eq. ([72]) for the class of initial states Eq. ([83]) . 



/•oo /"OO r 

p«(X,P) = dt dt7d3xdV^;(x',t',p)^«(x,t,p)Gi(t'-t)e'^**'(*-*)Qip(x,x',t-t'), (101) 
Jo Jo J 



where 



g^,p(x,x',.) = (<Ao|J-(x')Vnx,pe-VP^+^^(^)Jfl;,J+(x)|0o). (102) 



The amplitudes ^q,(x, i,p) in Eq. pU2|) are identical to the ones of Eq. ([M|) . however we have 
inserted explicitly their dependence on the mean momentum p of the initial state Eq. (j83p . We have 
assumed that the differences in initial momentum, between different mass eigenstates, are many orders 
of magnitude smaller than the sampling accuracy Sp, so that without loss of generality they may be 
considered equal. 

In Sec. 5.2.1 we showed that for position samplings the kernel Q^ p(x, x', s) decays to zero for values 
of s much larger than a suppression timescale Tsup- This time-scale arose from the diffusion of a state 
localized in an area of size 5, under the action of the free particle Hamiltonian. This time-scale is not 
affected by the momentum samplings because momentum commutes with the free particle Hamiltonian. 

We now consider that the detector lies at a distance L from the source, such that the scale 6tij for the 
time-of-arrival interferences is many order of magnitudes smaller than Ts^p. Hence, as in the derivation 
of Eq. ([93]) . we evaluate (5xp(x, x', s) at the limit s — )■ 0. We obtain 

Pa{^,P)= dt dt' d^xd^x'A*^{x',t',p)Aa{x,t,p)e''^^(^-^>q'^p{x,x'), (103) 

Jo Jo J 

where we wrote Q^p{x, x') = Qx p(x, x', 0). The kernel Q^ p defines a positive operator on the Hilbert 
space TiA of the oscillating particle. This operator does not correspond to a POVM because it is 
not normalized to unity. However, it does inherit the localization properties of the POVM IIxp- In 
particular, if |x — X| >> 5 and |x — X| >> 6, then Q^p{x,x') ~ since the scale 6 is macroscopic 
and hence much larger than the time scale associated to the microscopic currents. Moreover, from Eq. 
()102p . we see that if the detected particle's wave-function (as described by J'^ {x\cj}o))) is not localized 
on values of momentum within an accuracy Sp from P, the kernel Q" again vanishes. If the momenta 
of all product particles are measured, then the momentum p of the incoming particle can be inferred, 
by conservation. Hence, the action of Q" on the vectors of TIa leads to momentum localization with 
spread sp. If some of the product particles are undetected then the measurement does not allow the 
inference of the A-particle's momentum with the accuracy Sp, but with a larger spread s' depending on 
the details of the interaction. 

Consequently, Q" suppresses the values of the integral Eq. (jl03p . unless x and x' lie within distance 
of order 6 from X, and p lies within distance of order Sp from P. For the amplitudes in Eq. ()103p 
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there is no coherence for the position variable at the scale 6 (i.e., a « 5) and no coherence of the 
momentum variable at the scale of Sp (i.e., a~^ « Sp). Since Sp6 » h, this implies that Eq. (jl03p is 
well approximated by an expression of the form 

p,(X, P) = K{X, P) J d^x J d^pxsil^ - X|)x;,(|p - P|) (104) 

/ dt dt' / d3xd\U;(x,t',p)^„(x,t,p)e''*h(t-t)^ 
Jo Jo J 



X 



where x<5(^))Xs {p) ^-re positive functions localized around x = with width 5, and p = with width 
Sp, respectively. The multiplicative constant i^(X,P) has been inserted so that we can choose X5{x)^ 
Xg {p) to be normalized as 

Id^xx5ix) = l Jd^px5ip) = l. (105) 

Eq. (jlOSp shows that the effect of the uncertainties in the position and momentum samplings essen- 
tially correspond to a smearing of the probability distribution Eq. (|93l) , at a scale of 6 for position and 
Sp for momentum. If the relative sampling errors 6/L and (JpjP are much smaller than unity, then the 
oscillatory behavior of the probability density Eq. (I93p remains unchanged. 

5.4.2 General initial state 

In our previous derivation of the probabilities for particle oscillations we assumed a superposition of 
Gaussian initial states Eq. (I83p . However, in the Appendix B, we show that the initial state for 
particle oscillation experiments is a mixed state of the form Eq. (jl39p . that can be expressed as a 
convex combination of states p^oj C corresponds to quantities that are distinguished macroscopically in 
the production process. The macroscopic variables ^ allow for a coarse determination of the particles' 
position and momentum. We assume that the initial states of the target particles in the production 
region (see Appendix B) are quasi-classical so that no spatial quantum coherence is preserved in p^g- 
This means that we can express /5^o ^^ ^ combination of localized pure states such as Eq. ()83p . i.e., 

h^ = Xl ^fii^h / dqdpf^{q,p)\qp)ij{qp\, (106) 

where 

{x\qph = (^^^""^^"^ (107) 

the index i labels the mass eigenspaces. The function f^{q,p) is non-negative and it corresponds to a 
distribution of inferred values of position q and momentum p for the oscillating particle [119| . In Eq. 
(jl06p we have implicitly assumed that the uncertainties in the initial momentum p due to incoherent 
mixing are much larger than the differences \pi — pj \ of momenta in Eq. ()83p . 

Since the probabilities Eq. (jSip are convex functions of the initial states, we can write 

Po.iL) = f dqdpf^{q,p)p^^'P\L), (108) 



where pa'^ (L) are the probabilities Eq. (j93p computed for an initial state characterized by mean position 
and momentum q and p respectively. The Hamiltonian for free particles is translationally invariant, 
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hence, the variation of q can be absorbed in a redefinition of L — t- L + (7. For the regime where Eq. ([93 
holds, 



Pa{L) = / dqdpf^{q,p)K'{q,p) 



/ . \'-^l3i'-'ai\ e 



+ 5:2i?e('[/^,t/..C/^,[/:,.e-^^+^^^^+'^^e^^-(^+^)^ 



«<j 



(109) 



where fcjj is computed from Eq. (j90p for pj = pj = p. 

We note that the maximum scale of variation for q is set by the dimension d of the production region 
that satisfies d « L. We also assume that the uncertainty in the specification of the initial momentum 
Ap « p, where p is the mean momentum corresponding to the probability distribution f^{q,p)- The 
oscillating phases (pij = kij(L + g) is a variable of q and p. The mean value (p = {kij{L + q)) with 
respect to /^ approximately equals kijL, where kij is the oscillation wave-number evaluated for the mean 
momentum p. The corresponding spread in the phase is 

A4>ij < LAkij + Aqkj < 4>tj(.Y + —) « 4>ij- (HO) 

L p 

This implies that the substitution of the mean value of e**^'^ with e*''^'^ is a good approximatioro. Hence 
there is no significant modification in the oscillation phase when the effects of incoherent mixing is taken 
into account. 

To summarize, we have shown that Eq. (|90p for the oscillation wave-number remains valid for a 
large class of initial states of the form Eq. ()106p . The latter exhibit no spatial quantum coherence at 
scales larger than the width a of the Gaussian states Eq. (|83p . Hence, the effect of an initial state 
Eq. (1106P to the detection probability is formally similar to that of sampling uncertainties sampling, 
i.e., they both induce a smearing of the distribution probability by classical uncertainties. Thus, both 
sampling uncertainties and uncertainties from the initial condition contribute to the random errors of 
the measurement: they can only degrade the oscillation pattern of Eq. ()93p . and not modify its main 
features. 

5.4.3 Dispersion effects 

In order to examine the effect of dispersion on the probabilities, we calculate the amplitude Eq. ([83]) . 
by expanding Ei{p) to second order and Ti{p) to first order around pi, i.e., 

Ei{p) = Ei + Vi{p - Pi) + -fii{p - PifTiip) = Ti + Ai{p - Pi) (111) 

where 



Ei = \jnij+pl T, = T,{pi), Vi = ^, fi^ = ml/El A, = {dT , / dp){pi) . (112) 

We then obtain 

i^_M__ i{p,x-E,t)-V,t 



Ao,{t,x) =J2UpiUa 



\Ja'^ + iiXit 



^ The 'mean value' and 'deviation' of quantities such as (\>ij refers to the distribution of the parameters of Eq. (|83p in 
the weighting function. These objects describe features of the initial state and not statistics of measurement outcomes. 
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X exp 



"' ({x-v4-^? + ^t' + a^l^W 



2(^4 + ^2^2) ^ 

(113) 



+ 2(a4 + ;.2,2)l(^ -^* ^^ ) ^2 +A^* 



The effects of the dispersion become significant for time-scales such that a^ ~ ^jt. Since the ampU- 
tudes Eq. (|113|) are peaked for times t ~ L/vi the no-dispersion regime corresponds to the condition 

3^2 — << 1. The spread of the Gaussian o" is a free parameter that may take any value between the 

atomic and the weak timescale. For a ~ 10~^^m, Ei ~ lOOMey, the non-dispersive regime provides an 
adequate description if L « l^^^km. A value of a ~ lO^^^m is required, in order to see a failure of the 
dispersive regime at a distance L relevant to neutrino oscillation experiments. 

5.4.4 Particle detection through several processes 

We finally consider the case of a set-up where oscillating particles are detected by different processes. 
We show that the probability distributions for each process are added consistently into the expression 
for the total oscillation probability. 

Let us label by the index n = 1, . . . ,N the different channels through which oscillating particles can 
be detected. Each channel corresponds to a different Hilbert space T-Ln, that describes the non-oscillating 
particles participating in the n-th process. The construction in Sec. 5.2.1. then can be used with the 
following changes. 

1. The Hilbert space Tir is the tensor product ^n'Hn- 

2. The initial state \<po) of particles on the detector is of the form \(/)q ) ® . . . ® \(/)q ), in terms of 
states \(f>Q ) G Tin- 

3. The Hamiltonian Hr on the Hilbert space Tir splits as 

Hr = -/?(!) 0lO...(g)l-h... + lOlO...(g) i^(„), (114) 

in terms of Hamiltonian operators -ff(n) on Tin- 

4. The currents J^ (x) in Eq. (j67p for the interaction Hamiltonian are of the form 

J^(x) = J^^)^(x) 1 ® . . . ® 1 + 1 «) jP^(x) 0...01 + ... + 1^10...0 4^)^(x), (115) 

in terms of current operators Ja (x) on Tin. 

5. The positive operators H;^, corresponding to measurements of position of product particles, are of 
the form 

Ha = H^^ (» 1 . . . 1 + . . . + 1 ® 1 . . . H^f \ (116) 

in terms of positive operators IIx on Tin such that Hj^ |(/)(n)o) = 0. 
Then Eq. (|72p follows, with the kernel Rf{x,x',t—t') expressed as 



RU^,x',t-t') = E(4"^l4")-(xO\/n?e-^%)(*'-Vn?^i"'^+(x)l4"^ (117) 
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The remaining analysis of Sec. 5.2 passes unchanged, leading to 

p.{L) = E ^n ( E m^ar^'e-'^'' + E ^^^ (u^.U^.U^.K^e'^'^'-'^^'' e'''7' ^)\ ) . (118) 

n y i i<cj \ I I 

This is a convex combination of the probability densities Eq. (193p . where each process is weighted by 
different multiplicative coefficients K'^. The oscillation pattern of each term is modulated by wave- 
numbers k(^ as in Eq. (j90p . differing only by the values ej^ of the process' threshold. In other words, 
the detection channels are statistically independent and the resulting probabilities can be added. 

A Coarse-graining in time-of-arrival measurements 

In order to demonstrate how the decoherence time scale used in the definition of the POVM Eq. (j55p 
arises, we employ the formalism of Sec. 5.1 for a particle detection by a macroscopic apparatus. We also 
examine special cases for the POVM Eq. (j55p . showing that it recovers von Neumann measurements in 
a specific regime. 

We assume that the particle is annihilated or is absorbed by the apparatus, as part of the detection 
process. It is therefore convenient to employ the Fock space 

Ts = c®Us®{ns®ns)s®--- (119) 

associated to the Hilbert space Us of a single particle; S denotes symmetrization. We denote the Hilbert 
space describing the apparatus's degrees of freedom as Happ- The Hilbert space of the combined system 
is then Ts ® Tiapp- 

The Hamiltonian of the combined system is Hs(^l + l(^Happ + Hint, where Hs describes the dynamics 
of the microscopic system and Happ describes the dynamics of the apparatus. The general form of an 
interaction Hamiltonian describing the absorption of the microscopic particle in the process of detection 
is 

Hint = J2jd^x{ai{^)J+{^) + a|(x)jr(x)), (120) 

i 

where di(x),aj(x) are the annihilation and creation operators on J-g, i labels non-translational degrees 
of freedom (e.g., spin), and J^(x) are current operators on Tiapp with support in the region D where 
the detector is located. We assume that the microscopic system is prepared in the single-particle state 
l^o) and that the initial state of the apparatus |^o) is stationary. We set the scale of energy so that 
Happl'^o) = 0. Since the particle is absorbed by the detector we require that J~(x)|^o) = 0. 

The pointer variables correspond to a POVM 1 (8) 11;^, where Ha is a POVM on Tiapp defined so that 
/ dXIi\ = P is the projector on the excited states of the detector. In order to ensure that no detection 
events occur unless a particle enters the detector, we must assume that [Happ, P] = 0. It follows that 
[P,Hs'S>l + l(^ Happ] = 0. Hence, Eq. ^ apphes, 

C'A,t(l^o®)|^o)) = E/^'^(^^(^'*)|0))® (e'^"""*V^'^^WI^o)) , (121) 

where '(/'(x, i) is the evolution of iV^o) under e~*^°* in the position representation. 
Then, 



Tr [Cx,sPoCls') = E/'^'^^'^V^(^'^)V'K^''^')(^o|i-(x')\/ftIe-^^-''(^'-^)V^J+(x)|^&o). (122 
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In order to obtain a probability density p{X, t) we employ a smearing function \J fris — t) on both 
arguments s,s' of Eq. (jl22p : r corresponds to the temporal resolution of the apparatus. We note that 
the s and s' dependence in Eq. (jl22p arises from two independent terms: 'ipi{ii.,s) and the currents 
J"'"ib(x, s). In the former, the time dependence refers to the evolution of the microscopic system prior 
to detection, while in the latter it refers to properties of the macroscopic apparatus. The temporal 
resolution of the apparatus is determined by the latter terms. 

In order to identify the temporal resolution of the apparatus, we examine the expectation value 



P{s) = (^o|J~(x )VnAe-^^''-VnAi+(x)|vI'o) (123) 

from Eq. (jl22p . This satisfies 

\P{s)\ < K^o|i7(x )\/ftIe-*^"-^V^J+(x)|vI/o)| • ||<%-*^--^<'||2, (124) 

where || • II2 stands for the Hilbert-Schmidt norm. The s-dependent terms in the right-hand side of Eq. 
()124|) equal jTr ( Jllxe~^^'^pp^ylixe^^'^pp^ ) and they can be written as 



i^l^ f!l ^ (125) 

\ TrUx 



in terms of a density matrix p = ylix/ \lTrIi\. 

The object Tr[e "pp^ pe~'^ "pp^ p] is known as the persistence probability of the state p; its properties 
have been studied in relation to exponential decay J114j . Typically, a persistence probability is charac- 
terized by a decay time-scale Tdec, for s >> Tdec the persistence probability becomes negligible. In the 
context of Eq. (|122p . this implies that interference between terms with different values of s and s' in 
Eq. (I122|) are suppressed if |s — s'| >> Tdec- In particular, this implies that the interference terms in 
Eq. (I50p become negligible if the width of the intervals is much larger than Tdec- Restricted to such 
intervals, Eq. (|49p defines a probability measure. Equivalently, one may employ Eq. (j54p with functions 
fr characterized by values of the smearing parameter r >> Tdec- 

In order to estimate the decay time-scale Tdec, one should consider specific models for the self- 
dynamics of the pointer variable. In quantum Brownian motion models, for example, one describes the 
pointer variable as the position of an harmonic oscillator, of effective mass ;U*, in contact with a heat 
bath at temperature T. Taking the positive operators IIa to correspond to a sampling of the pointer 
variable with an accuracy of order 6, we find that Tdec ~ Trei{fJ'*T6'^)~^ , where r^ei is the relaxation time 
for the pointer variable [43]. 

In conclusion, a choice of r in the smearing functions of the order of Tdec, ^-s determined by the 
dynamics of the apparatus guarantees that the interference terms in Eq. (fSOJ) are negligible. Hence, Tdec 
sets a scale of minimal temporal coarse-graining necessary for the quasi-classical attribution of definite 
macroscopic properties to the pointer variables. 

Special cases. If the states ipi{x, s) depend on s at a scale much smaller than the characteristic time- 
scales in J tl\e~'^^'^pp^ y Ii\, then we can approximate 'Jti\e~^^'^^ ^^W 11^ ^ Ha in Eq. (|122p . In this 
regime, the time of arrival cannot be distinguished. Therefore, we consider the partial probabilities p{\) 
obtained by integrating p{t,X) in Eq. (|54p . over the time of arrival. Then, Eq. (j6Up applies and we 
obtain 

pW = ^E/^'^^'^'^^j(^)(^'^') (l^dsi;i{x,s)] (£ds'x(jj{x',s')\ , (126) 



V 
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where Q(A) is a family of positive operators on the Hilbert space T-Ls with matrix elements 

(x'|Qi,-(A)|x) = 27rr2(xI/o|J-(x')nAJ+(x)|^o>, (127) 

where the factor 27rr^ was inserted into the definition of Q so that {ijj\Q{X)\ijj) has the correct dimensions 
for a probability density over A. 

If, on the other hand, time-scales in "01 (x, s) are much larger than Tdec, we can substitute the smearing 
functions \/J{t) in Eq. (j52]) with delta- functions, i.e., we set yj frit — s) = (27rr^)^/^5(t — s). We then 
obtain 

p{t,\) = -jL=Y, j d^xd^x'^i{^,t)i)j{-s!,t)Qi,{\){-K,^). (128) 

The function p[t, A) in Eq. (jl28p is a proper density with respect to the time of arrival t, so it can 
be integrated over t yielding a probability density p{\) 

^(^) = -7^ r dt{i^t\Q{X)\i^t). (129) 

V27rr^ JO 

where |V't) = e'^-'^^'I'^o)- 

In Eq. ()129p the probability p(A) is defined as an integral over the single-time probabilities ('(/'i|Q(A)|'0i) 
As we also argued in Sec. 4.2.3, this result applies only to the specific regime we considered here. 

Von Neumann measurements can be recovered from Eq. (|128p . if the initial states are sufficiently 
localized. This means that the distribution of the time of arrival should be sharply peaked around 
classical arrival time tcu with a spread indistinguishable by the temporal resolution of the detector. 
Indeed, setting t = t^i in Eq. ()128p . we obtain p{X) ~ ('(/'t^JQ(A)lV'tci)- 

When we integrate p{t,X) over A in Eq. (jl28p . we construct a probability density for the time of 
arrival t, 

p{t) ~ {i^timt), (130) 

where B = f d\Q{\). In this regime, 13 acts as a probability current operator. Its explicit form depends 
on the interaction between the microscopic system and the detector, i.e., on the currents J±. We also 
note that in the regime under consideration, decoherence has rendered the time of arrival a quasi-classical 
variable, so that the expectation value of a current operator is positive- valued — see Refs. [56] . 
We can approximate B with an operator whose Wigner function is of the form 

W{x,p) = f{\p,\)6{x-L), (131) 

by assuming that the size d of the detector is much smaller than the source-detector distance L and that 
the detector is homogeneous, so that for x and x' lying in the detection region, {'i>Q\J-'{x.')J^{x.)\'ip()) is 
a function of |x — x'|. 

The positive function / in Eq. (I13ip depends on the details of the interaction. The standard 
probability current (modulo terms of higher order in h) corresponds to f{\px\) ~ \Px\- For / a constant, 
we obtain a particularly simple form, p{t) ~ \ilj{L,t)\'^. 
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B The role of the initial state in neutrino production 

In this section, we re-derive a fundamental property of quantum measurement theory, in a context 
relevant for neutrino oscillations. We consider the most general set-up for neutrino (or any particle) 
production and we show that all physical information that can be extracted by measurements is encoded 
in the reduced density matrix of the neutrinos. One reason for including this derivation is that a large 
number of articles associate the standard expression of oscillation wavelength with special initial states 
of the neutrinos. The results in this section allow us to argue that such special conditions are highly 
implausible in a treatment of the production process. Another reason is that this derivation demonstrates 
conclusively why a probability assignment, linear with respect to the initial density matrix, is necessary 
in any precise quantum mechanical treatment of particle oscillations. 
We consider a production reaction of the form 

a + b^iy + di + d2 + ..., (132) 

where a are the particles of an incoming beam, b are the target particles, and di product particles of 
the reaction other than the neutrinos. There are many 6-type particles in the target, but each neutrino 
is produced in a reaction that only one of them participates. Hence, to each production event we can 
associate a unitary operator 

V -.Tia'^nb^nu'^nd.'^nd^'^..., (133) 

where 7ia,Tib, ■ ■ ■ are single-particle Hilbert spaces for the a,b, . . . particles. 

We define the corresponding Fock spaces Ta = C (B Tia ® (T-la '^ T~ia)s,A © , and similarly for 

other particles; S ox A stands for symmetrization or anti-symmetrization respectively. The total Hilbert 
space of the system is T-Ltot = ^a® ^b® ^u ® ^di ® ^d-x ® ■ ■ ■■ 

We typically assume that the reaction of a single particle in the beam corresponds to a single element 
of the statistical ensemble described by the quantum state. Hence, the a-particles are initially found in 
a single-particle state |V'a)a- The 6-particles are assumed initially in an A^-particle state |0i, . . . 4>N)h- 
The initial state for the neutrinos and the di particles is the vacuum. Hence, the initial state for the 
system is 

IV'a) © Ic/^i, . . . , 07V) © |0)^ ® |0)d, ® |0)d2 © . . . (134) 

An a-particle interacts through the channel Eq. ()133p with one of the b particles, thus the quantum 
state after the interaction is 

1^) = |0)a ® E (I'^l' •••,#.>••-, ^n)i. ® [V(|^)a ® l0.)6)]Mi,d2,...) ' (^^S) 

i 

where |0i, . . . ,^j, . . . 0Af)f, is the {N — l)-particle state obtained by removing the i-th entry from the 
A^-particle state | (/>!,... , (/>Ar). 

Reactions of the form Eq. (|132p can often be identified macroscopically by the detection of some 
of the produced particles dj, and by the ensuing measurement of some of their properties (position, 
momentum, and so on) collectively denoted by ^. The variables ^ may be defined at different moments 
of time and they may include the times of detection for the di particles. In Ref. [120] . it is argued that 
the neutrino and the other produced particles are initially entangled and that a detection of the latter, 
causing disentanglement, is a necessary condition for neutrino oscillations. 

Let n^ be a POVM on Utot that corresponds to measured properties of the di particles, as above. 
Then the ensemble of all neutrino production events splits into sub-ensembles characterized by different 
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values of ^. Each sub-ensemble is described by the density matrix 



= V ^' :^ 'V ^ (136) 

Subsequent evolution of p^ proceeds under a unitary operator JJtot on litot for the dynamics of 
all degrees of freedom. To leading approximation, there is no interaction between neutrinos so JJtot 
factorizes as Uy ®Ur, where Ur is defined on the Hilbert space of the remaining degrees of freedom 
T~ir = J^a ® J^b ® ^d-i ® ^d2 ® ■ ■ ■■ Then, for any subsequent measurement on the neutrinos, the full 
information for this sub-ensembles is contained in the reduced neutrino density matrix 



Tr-Hr 



U^ Ur)Ju^\^){'^\Ju^iUi U}) 



Tr-Hr 



iiA^){^dn 



(*|n$l*> (^in^l*) 



■Ul (137) 



The probabilities for the full ensemble will be obtained from the convex combination of p^^^ weighted 
by the probabilities (^|n^|^) corresponding to each sub-ensemble. Hence, the probabilities for any 
measurements of the neutrino will be determined by the reduced density matrix 

Pu = UuPuoUl (138) 

where 

p,o = J2'^rn^ (yri^l^X^iyS^) (139) 

is the initial reduced density matrix of the produced neutrinos. 

We showed therefore that, whatever the neutrino production may be, any subsequent measurement 
on the neutrinos depends only on the reduced density matrix of the neutrinos that evolves under the 
neutrino self-dynamics. In other words, the details of the particle production are all encoded in the 
initial density matrix, the same way that the details of the detection are encoded in a POVM. Moreover, 
the existence of any macroscopically retrievable information from the production process, implies the 
existence of different macroscopically distinguishable sub-ensembles. Thus, p^ is necessarily a mixed 
state. 

We also note the importance of having a probability assignment linear with respect to the initial 
state pi^: without linearity it is impossible to obtain a consistent combination of different macroscopically 
distinguishable sub-ensembles. This means that any probability assignment must be linear with respect 
to the initial state, otherwise it cannot provide the basis of sound probabilistic reasoning. This point is 
important in the critique of existing approaches to particle oscillations. 
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